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EVALUATION 


This  effort  provides  an  in-house  capability  to  demonstrate  advanced 
techniques  in  off-line  orthophoto  production.  Exploitation  of  multiple 
digital  terrain  data  sources  has  the  potential  to  greatly  expand  the  Defense 
Mapping  Agency  orthophoto  production  capability.  Orthophotos  are  a viable 
form  for  photogrammetric  targeting  data  bases  which  are  currently  being 
supported  by  TPO  2.  Positional  data  bases  in  an  orthophoto  format  require 
minimal  field  exploitation  equipment  complexity.  Orthophotos,  however,  are 
relatively  expensive  to  produce  due  to  the  high  cost  of  digital  terrain 
compilation.  By  taking  advantage  of  terrain  data  already  in  existance,  the 
need  for  new  compilation  with  new  image  coverage  is  negated.  While  contrac- 
tual resources  for  complete  testing  were  not  available,  this  effort  has  pro- 
vided experimental  software  which  will  allow  future  evaluation  of  this  potential. 

Project  Engineer 
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1. 


INTRODUCTION 


This  report,  submitted  in  fulfillment  of  RADC  Contract  No. 
F30602-74-C-0297,  presents  a description  of  the  PDP-11/45  software 
developed  by  DBA  Systems,  Inc.  for  the  offline  processing  of  digital 
relief  data  on  the  Offline  Orthophoto  printer  (OOP).  The  primary 
objectives  of  this  contract  consisted  of  the  following  items: 

a)  to  eliminate  the  dependence  of  the  OOP  on  the  AS-llB-1, 

b)  to  make  use  of  available  digital  data  base  information 
on  the  OOP, 

c)  to  allow  for  more  flexible  use  of  data  created  by  the 
AS-llB-1, 

d)  to  allow  for  creation  of  orthophotographs  from  limited 
relief  information  (discrete  random  points). 

The  Offline  Orthophoto  P^rinter  Software  (OOPS)  contains  a new 
approach  to  offline  processing  of  digital  data.  The  OOFS  accepts  in- 
puts consisting  of  terrain  data  and  photographic  orientation,  performs 
a series  of  computations,  and  outputs  two  data  files.  These  two  files 
duplicate  the  format  of  those  generated  by  the  AS-llB-1  stercoplotter 
and  can  be  directly  input  into  the  Offline  Orthophoto  printer  (OOP). 

The  first  file,  the  ^hutDown  Tape  (SDT),  contains  the  information 
necessary  for  orientation  of  the  photographic  imagery  on  the  OOP.  The 
second  file,  the  OOPS  Profile  Tape,  contains  the  derived  terrain  infor- 
mation necessary  to  produce  a rectified  orthophotograph. 


A detailed  description  of  the  software  written  under  the  OOPS 
contract  is  presented  in  the  next  eight  sections.  Section  two  contains 
a system  description.  This  includes  an  overview  of  the  previous  means 
of  generating  an  orthophoto  as  well  as  a description  of  the  new  approach 
developed.  Section  three  presents  a detailed  description  of  the  methods 
implemented  to  process  available  digital  topographic  data.  This  des- 
cription includes  profile,  contour  and  random  point  data  processing 
along  with  the  processing  of  multiple  sets  of  data.  Section  four  con- 
tains the  mathematical  concepts  used  to  develop  the  OOPS.  Section  five 
contains  the  user's  input  description.  Section  six  presents  the  system 
and  subroutine  flew  diagrams.  Section  seven  contains  a short  description 
of  all  subroutines  used  in  the  software  development. 

The  input  formats  of  the  terrain  data  sources  is  addressed  in 
Section  eight.  A description  of  the  two  output  files  is  contained  in 
Section  nine.  Results  of  various  data  runs  are  described  in  Section 
ten.  A summary  containing  conclusions  and  recommendations  is  presented 
in  Section  eleven. 
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2. 


SYSTEM  DESCRIPTION 


The  OOPS  system  has  substantially  increased  the  production 
capability  of  the  orthophoto  printer  at  RADC.  Before  the  implementation 
of  this  software  system,  the  OOP  was  totally  dependent  on  the  AS-llB-1 
stereoplotter  for  the  generation  of  all  input  data.  OOPS  was  designed 
to  expand  the  capability  of  the  OOP  by  accepting  various  types  of  digital 
data  for  the  computation  and  output  of  profiles  in  AS-llB-1  format.  The 
software  includes  options  to  input  orientation  data  in  various  coordinate 
systems,  transform  such  data  into  the  appropriate  coordinate  system  for 
generating  a punched  shutdown  tape  for  input  to  the  OOP. 

The  OOPS  system  is  written  primarily  in  Fortran  IV  but  includes 
some  Assembly  language  routines  for  special  I/O  handling.  The  system 
operates  on  RADC's  PDP  11/45  in  the  32K  of  available  core  using  the  DOS 
operating  system.  Available  peripherals  include  a card  reader,  paper  tape 
reader/punch,  two  seven  (7)  track  and  one  nine  (9)  track  magnetic  tape 
drives,  and  an  80  column/line  printer. 

2. 1 System  Overview 

The  capability  of  generating  an  orthophoto  has  significantly 
changed  with  the  creation  of  the  OOPS  system.  It  is  no  longer  necessary 
for  the  input  photograph  to  be  one  of  the  stereo-pair  used  to  generate  the 
digital  data  tape.  The  terrain  information  can  now  be  in  the  form  of 
cartographic  data,  AS-llB-1  profile  data,  or  discrete  points,  while  the 
photographic  orientation  parameters  can  be  derived  from  block  or  strip 
triangulations,  resections  or  AS-11  absolute  model  data. 


Previously  the  OOP  was  entirely  dependent  on  the  AS-llB-1 
stereo- pi  otter  for  the  generation  of  orthophotos.  Figure  2.1  shows 
the  means  by  which  an  orthophoto  was  produced  on  the  OOP  before  the 
implementation  of  the  offline  orthophoto  software.  This  figure  clearly 
demonstrates  the  dependence  on  the  AS-llB-1  for  all  data  input  to  the 
OOP.  In  addition  to  dependence  on  input  data,  the  input  photograph  was 
required  to  be  one  of  the  stereo-pair  used  to  generate  these  data  tapes. 

Figure  2.2  represents  an  overview  of  the  OOPS  system.  Com- 
parison of  Figure  2.2  to  Figure  2.1  shows  the  flexibility  of  input  data 
resulting  from  the  OOPS.  The  OOPS  system  will  accept  terrain  data  for- 
matted as  MMS-32  Cartographic  data,  AS-llB-1  profile  data,  or  discrete 
points,  such  as,  triangulation  data  and  survey  or  geodetic  control. 
Combinations  and  multiple  sets  of  the  above  terrain  data  are  also  ac- 
cepted. The  source  imagery  accepted  by  OOPS  includes  both  panoramic  and 
frame  photography.  Position  and  orientation  data  generated  for  this 
source  material  can  be  input  in  nearly  any  form.  Optional  coordinate 
systems  include  USR,  LSR,  LV,  geographies,  and  AS-11  absolute  model  para- 
meters. 

2.2  System  Computational  Flow 

Figure  2.3  presents  the  four  major  sequential  elements  used  in 
the  computation  of  simulated  AS-llB-1  data  necessary  for  orthophoto  gen- 
eration. These  four  elements  are  input,  extraction  of  data,  computation 
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FIGURE  2.1,  Previous  Means  of  Generating  an  Orthophoto  Printer  Product 


SOURCE  IMAGERY 


Extract  Area  Sort  Data  Computations 


FIGURE  2.3.  OOPS  (Off-line  Orthophoto  System) 


of  functional  terrain  representation  and  output  of  newly  created  pro- 
files to  a tape  formatted  similar  to  that  created  on  the  AS-llB-1. 

The  input  element  contains  the  logic  to  control  the  input  data  necessary 
for  the  computation  of  orientation  and  selection  of  the  correct  terrain 
data  link.  An  explanation  of  the  many  input  options  follows  in  Section 
4.1. 


The  general  concept  for  extracting  data  fromthe  source  tape  is 
the  same  in  all  modes.  Basically,  the  geographical  area  of  interest  is 
determined  by  projecting  the  film  format  corners  to  the  ground  and  then 
transforming  these  ground  coordinates  to  digital  data  coordinates.  This 
computation  defines  the  area  from  which  data  is  to  be  extracted.  As  the 
terrain  data  tape  is  read,  only  data  within  the  desired  area  is  stored 
on  random  disc,  thus  efficiently  using  the  random  storage  routine.  When- 
ever possible,  data  sorting  is  accomplished  during  this  storage  process 
to  eliminate  unnecessary  I/O.  The  sorting  algorithms  for  each  data  type 
are  explained  in  detail  in  Section  3. 

Computation  of  the  relief  model  and  photographic  orientation  are 
two  of  the  most  important  areas  in  the  OOPS  program.  Although  each  data 
type  is  handled  differently  due  to  considerations  of  program  efficiency 
and  accuracy,  the  basic  computational  objective  is  the  derivation  of  a 
mathematical  representation  of  the  terrain  surface.  A description  of 
these  mathematical  functions  is  presented  in  Section  3.  Other  computa- 
tions required  include  transformation  of  the  input  photo  orientation 
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and  position  parameters  to  the  coordinate  system  of  the  OOP.  These 
transformations  are  explained  in  Section  4.1. 

The  last  major  area  of  OOPS  deals  with  data  output.  The  shut- 
down tape  contains  the  necessary  information  for  photo  orientation  along 
with  values  for  atmospheric  refraction,  earth  curvature,  focal  length 
and  various  other  values  needed  on  the  BX-272.  The  shutdown  tape  is 
punched  on  paper  tape  in  RCA  501  code.  Hard  copy  output  of  the  shut- 
down tape  information  is  tabulated  in  a convenient  format. 

The  only  other  major  output  of  the  OOPS  program  is  the  simu- 
lated AS-llB-1  profile  tape.  This  output  tape  is  formatted  as  "X,Y,Z" 
coordinates  in  the  LSR  system  required  by  the  OOP  system.  The  detailed 
format  of  this  tape  is  documented  in  Section  8. 
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3. 


METHODS  OF  PROCESSING  VARIOUS  DATA  TYPES 


The  OOPS  system  will  accept  various  types  of  digital  relief 
information  and  derive  profiles  in  the  desired  model  coordinate  system. 
The  methodology  selected  for  processing  a specific  data  type  resulted 
from  considerations  of  program  efficiency  and  resultant  accuracy.  The 
following  sub-sections  describe  in  detail  the  method  of  processing 
each  data  type. 


An  algorithm  based  upon  the  multi  quadric  function  is  used  to 
model  the  terrain  relief  in  the  profile  and  random  point  modes,  while 
a polynomial  function  is  used  in  the  contour  mode.  These  algorithms 
will  be  presented  at  this  level  in  order  to  clarify  subsequent  sections. 
The  basic  equation  for  the  multiquadric  function  is  of  the  form  devel- 
oped by  Hardy  (see  Reference  1)  and  is  given  by: 


where 


1. 

z 


Z = L a,/((Xi-X,)s+(Y, -Yj)2+C)  . 
1=1 


(3.0.1) 


Xi,Y, 

C 


Z coordinate  of  jth  data  point  in  local  space 
rectangular  system 

coefficient  associated  with  the  ith  nodal  point 

coordinates  of  ith  nodal  point 

coordinates  of  jth  data  point 

an  arbitrary  constant  term  representing  the  area 
of  influence  of  nodal  points. 
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It  is  the  unique  manner  in  which  Hardy's  formulation  is  applied  that 
is  important  in  the  OOPS  system.  This  application  is  described  in 
Section  3.1  and  3.3. 

The  equation  for  the  polynomial  representation  of  relief  in 
the  contour  mode  is  given  by: 

Z = + QgX  + 03/+ + ag/S  + (3.0.2) 

where 

X = X - X 

y = y - y 

n 

X = L Xj 

1 = 1 

y = S y, 

1 = 1 

n = tot-al  number  of  data  points  used  in  a given  raster  fit 
x,y  = model  coordinates. 

The  number  of  polynomial  coefficients  used  in  a specific  raster 
fit  is  dependent  upon  the  number  of  available  data  points. 

3. 1 Profi le 

Data  input  to  the  profile  mode  is  currently  derived  from  the 

AS- 11  plotter  system.  This  data  is  formatted  in  a very  systematic 

manner.  The  data  is  very  dense  and  is  formatted  as  regularly  spaced 

"X,Y"  model  coordinates  with  associated  "Z"  values  in  a local  earth 
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tangent  coordinate  system.  Advantage  is  taken  of  the  systematic  for- 
mat of  this  data  in  a very  unique  manner. 


The  partial  derivative  of  equation  3.1  with  respect  to  the  un- 
known coefficients  is 

1 

6Z/60j  = I.0/(X,  -Xj)2  +(Y,  - + C)^  (3.1.1) 

Thus,  the  partial  derivative  is  simply  a function  of  the  distance  of 
the  jth  data  point  from  the  ith  nodal  point.  With  evenly  grided  data 
this  is  an  important  criteria,  for  if  one  selects  subsets  of  the  grided 
data  and  establishes  an  origin  at  the  center  of  each  subset,  the  data 
becomes  perfectly  symmetrical.  This  means  that  if  a consistent  nodel 
point  grid  is  designated  within  every  data  subset  the  partial  deriva- 
tives computed  using  equation  (3.1.1)  are  identical  for  the  jth  data 
point  in  all  subsets.  This  implies  that  the  least  squares  normal 
equations  for  all  subsets  are  identical  and  only  the  constant  vector 
changes  as  a result  of  the  changes  in  elevation  within  each  subset. 

As  a matter  of  fact,  one  can  form  and  invert  the  normal  equations  and 
reuse  the  same  inverted  normal  equations  to  solve  for  coefficients  in 
all  subsets  which  have  equivalent  spacing.  The  OOPS  system  currently 
employs  this  technique  in  the  profile  mode.  The  inverted  normal  equa- 
tions are  formed  and  stored  for  only  one  subarea  of  the  total  data  set 
and  reused  for  all  subsequent  subareas.  In  this  regard,  future  versions 
of  the  program  could  be  designed  to  accept  as  input  one  or  more  sets  of 
inverted  normal  equations. 
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Figure  3.1.1  shows  the  method  of  subarea  division  implemented. 

The  origin  of  each  subarea  is  designated  by  the  large  solid  triangles. 

The  wide  dashed  lines  surround  data  areas  that  are  stored  as  a unit. 

New  profiles  are  to  be  derived  over  thr  area  bounded  by  the  dotted  lines. 

Although  the  multiquad  solution  is  derived  over  four  of  the  squares  out- 
lined by  the  wide  dashed  lines  this  representation  is  only  used  over 
the  inner  subarea  outlined  by  the  solid  black  lines  with  the  solid 
black  triangles  as  the  center.  The  excess  area  is  included  in  the  least 
squares  solution  to  avoid  discontinuity  problems  between  adjacent  areas. 

i 

The  total  area  is  configured  as  a square  to  simplify  the  computational 
algorithms. 

In  the  profile  mode,  the  digital  data  is  handled  in  four  basic 
steps.  First  the  software  determines  the  limits  of  the  required  area; 
next  this  required  data  is  stored  in  a systematic  manner;  the  third  step 
entails  the  derivation  of  the  multiquadric  coefficients;  and  the  final 
step  is  the  derivation  of  the  new  profiles  using  the  derived  coefficients. 

Initially  the  program  software  derives  the  object  space  bL^nd- 
aries  of  the  relief  data  required  to  produce  the  new  orthophotograph 
(see  Section  4 for  coordinate  system  definition).  These  boundaries  are 
determined  by  projecting  the  corners  of  the  user  designated  photo  area 
to  the  earth's  ellipsoid.  These  projected  corners  are  referred  to  as 
the  footprint  of  the  required  orthophoto  area. 
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Origin  of  subareas 

Data  subsets  stored  as  units 

Area  new  profiles  are  to  be  derived  over 


Area  coefficients  derived  over 
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The  digital  terrain  data,  as  required  to  represent  the  projected 
footprint,  is  read  from  magnetic  tape.  The  required  data  is  stored  in  a 
systematic  manner  on  random  disc.  Figure  3.1.1  shows  the  subareas  of 
digital  terrain  data  which  are  stored  as  units. 

The  next  step  involved  in  the  generationof  profi le  data  is  the 
actual  formation  of  the  normal  equations  and  solution  for  the  multi- 
quadric coefficients  describing  the  terrain.  Nodes  are  computed  over 
an  area  that  is  two  times  the  size  of  the  coefficient  block  in  both 
the  X and  y direction.  These  nodes  are  evenly  spaced  about  the  center 
of  this  square  area.  The  square  area  is  four  times  the  size  of  the 
storage  block  to  allow  for  a fifty  percent  overlap  on  all  sides. 

The  normal  equations  are  then  formed.  The  equation  used  is  a 
function  of  the  x and  y coordinates  (see  equation  3.0.1).  Recall,  be- 
cause of  the  grided  format  of  the  profile  data,  the  least  squares  normal 
equations  are  the  same  for  all  subareas.  Therefore  the  inverse  for  the 
normal  equations  is  stored  on  disc  after  the  first  solution  and  used 
throughout  the  digital  data  area.  Only  the  constant  vector  is  recomputed 
for  subsequent  blocks.  Having  recomputed  the  constant  vector,  the  multi- 
quadric coefficients  describing  the  terrain  are  derived  (Section  4.3)  and 
stored  on  a disc  file  for  later  use  in  computing  the  final  output  pro- 
files. In  order  to  evaluate  the  adequacy  of  the  derived  terrain  function, 
a RMS  value  between  the  derived  "Z"  coordinates  and  the  input  profile 
values  is  computed.  This  is  computed  by  deriving  "Z"  values  using  the 
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multiquadric  function  at  all  data  points  of  every  profile  used  for 
that  particular  solution.  The  difference  between  the  original  value 
and  the  derived  value  is  used  to  compute  the  MEAN  and  RMS  for  each 
specific  subarea.  This  statistic  is  output  so  the  user  can  determine 
the  quality  of  the  fit  over  each  subarea.  Also  output,  are  summary 
statistics  which  derive  the  MEAN  and  RMS  values  for  the  entire  area. 

With  all  the  multiquadric  coefficients  computed  and  stored  on 
a disc  file,  the  last  step  in  deriving  the  new  profiles  begins.  The 
coefficients  for  a subarea  are  retrieved  and  used  to  derive  the  "Z" 
coordinates  for  "x,y"  values  within  the  bounds  of  that  subarea.  These 
derived  profiles  are  reformatted  to  correspond  to  the  AS-11  format  in 
the  output  link  (see  Section  9). 

3. 2 Contour 

One  of  the  objectives  of  the  OOPS  system  is  to  exploit  data 
from  the  Automatic  Cartographic  System  (ACS).  The  format  of  the  digital 
data  currently  used  in  the  ACS  has  been  designated  as  MMS-32.  The  ter- 
rain relief  data  used  in  the  ACS  is  stored  as  contour  data.  Thus,  the 
MMS-32  contour  data  does  not  have  the  symmetrical  "x,y"  grid  pattern 
of  the  profile  data.  It  does,  however,  have  an  important  characteristic. 
As  terrain  slope  increases,  a greater  number  of  contours  are  used  to 
represent  an  area.  That  is,  relief  information  is  denser  in  areas 
characterized  by  significant  changes  in  ground  elevation. 
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The  following  general  approach  for  deriving  profiles  from  the 
MMS-32  digital  contour  data  has  been  implemented.  The  final  model  sys- 
tem is  established  (see  Section  4)  and  the  segment  of  contour  data  near 
a desired  profile  is  stored  as  "raw"  information  for  that  profile.  The 
"raw"  profile  information  is  discrete  points  which  are  taken  from  the 
digital  contour  data.  These  extracted  discrete  points  are  those  nearest 
the  intersection  of  a contour  and  the  line  along  which  a profile  is  to 
be  derived.  Desired  profiles  are  those  which  are  to  be  created  over  a 
user  specified  photo  area.  These  profiles  have  a constant  "x"  coordin- 
ate in  the  newly  established  model  system.  The  "y"  coordinates  are  in- 
crements of  250pm,  while  the  "z"  coordinate  represents  the  terrain  ele- 
vation above  an  established  datum.  When  all  the  contour  information  has 
been  stored  on  random  disc  in  this  manner,  the  "raw  profiles"  are  pro- 
cessed in  sets  of  four.  Polynomial  functions  are  fit  over  raster  areas 
covering  the  width  of  four  profiles,  but  vary  in  length  depending  on 
the  density  of  the  data.  Thus  the  size  of  the  area  analytically  repre- 
sented is  directly  dependent  on  the  terrain  slope.  Following  is  a de- 
tailed description  of  the  processing  sequence  performed  by  the  contour 
link. 


The  MMS-32  tape,  which  contains  cultural,  hydrographic,  and 
topographic  data,  is  searched  until  a data  block  containing  contour 
data  is  found.  Each  data  point  within  that  block  is  then  transformed 
to  the  desired  profile  system  and  checked  to  see  if  it  lies  within  the 
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desired  model  limits  of  the  designated  orthophoto  area.  All  data  points 
which  are  within  the  required  limits  are  stored  along  with  their  profile 
identification.  The  available  core  storage  is  divided  into  blocks  ac- 
cording to  the  total  desired  number  of  output  profiles.  Thus,  a parti- 
cular set  of  profile  data  is  stored  in  a block  of  core  (see  Figure  3.2.1). 
For  example  profiles  1-5  may  be  stored  in  block  1 while  profiles  6-10 
would  be  stored  in  block  2,  etc.,  through  block  19.  Block  number  20 
contains  the  profiles  which  are  not  evenly  divisible  into  the  first  19 

blocks.  On*  proftle  set  stored  (n  each  block 
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FIGURE  3.2.1.  Data  Storage  (Contour) 
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As  each  block  fills,  it  is  transferred  to  random  disc.  Each  series  of 
blocks  are  chained  together  on  random  disc  by  including  the  address  of 
the  previously  written  block  in  the  subsequent  block.  Thus,  all  blocks 
containing  a particular  subset  of  profiles  can  be  retrieved  by  knowing 
the  address  of  the  final  block.  Each  subset  of  profiles  is  retrieved 
and  the  profiles  are  written  on  a sequential  file  in  order  of  ascending 
"y"  coordinates.  If  more  than  one  point  is  found  for  the  specified  "y" 
profile  interval,  only  the  point  closest  to  the  desired  profile  incre- 
ment is  retained. 

The  profiles  are  retrieved  in  blocks  of  four  and  polynomial  fits 
are  computed  over  small  raster  areas  of  the  four  consecutive  profiles. 
Using  a set  of  four  profiles  the  two  center  profiles  are  derived  and  out- 
put for  each  raster  area.  This  concept  is  illustrated  in  Figure  3.2.2. 
For  a normal  raster,  this  is  accomplished  by  deriving  a least  squares 
solution  for  the  seven  polynomial  coefficients  in  equation  3.0.2.  The 
data  points  used  for  this  solution  are  the  raw  profile  points  as  des- 
cribed earlier  in  this  section,  plus  the  two  end  points  derived  during 
the  previous  raster  fit.  A normal  raster  is  defined  as  one  not  contain- 
ing the  end  of  a profile  set.  Thus,  ten  (10)  new  data  points  are  avail- 
able and  constitute  a raster  area.  For  nonnormal  rasters,  less  than  ten 
(10)  new  data  points  are  available,  the  number  of  coefficients  used  in 
the  least  squares  fit  is  reduced.  If  four  (4)  through  nine  (9)  new 
points  are  available  five  (5)  coefficients  are  used,  if  one  (1)  through 
three  (3)  new  points  are  available  three  (3)  coefficients  are  used. 
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After  the  initial  fit  the  two  end  points  from  the  profiles  derived  in 
the  previous  raster  are  retained  and  included  as  data  in  the  subsequent 
raster  fit.  These  data  points  are  retained  to  assure  that  consecutive 
rasters  fit  together  properly.  Each  of  these  retained  data  points  are 
weighted  as  tl.ree  points  to  assure  raster  compatibility.  Consistency 
between  profile  sets  is  accomplished  by  saving  two  profiles  from  the 
previous  set  and  utilizing  only  two  new  profiles. 

These  algorithms,  implemented  in  the  OOPS  contour  link,  pro- 
vide an  efficient  and  accurate  method  of  deriving  profiles  from  the 
MMS-32  contour  data.  Using  this  set  of  four  profiles  two  output  pro- 
files are  derived  along  the  center  two.  This  concept  is  illustrated 
below. 


Raster 

Area 


Raster 

Area 


Derived  profile  coordinates 

A Digital  data  points  as  extracted 
from  MMS-32  contour  tape 

FIGURE  3.2.2  Raster  Area 
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The  above  method  for  deriving  profiles  using  contour  data  pro- 
vided excellent  results  over  areas  of  moderately  sloping  terrain.  How- 
ever, over  areas  of  steep  terrain,  the  polynomial  function  did  not  com- 
pletely conform  to  the  terrain.  To  more  accurately  model  the  steep 
terrain,  an  interpolation  routine  was  developed  from  which  an  elevation 
could  be  computed  based  on  the  distance  the  desired  data  point  was  from 
the  known  data  point.  The  distance  from  the  desired  data  point  to  the 
known  data  point  is  computed  from: 


dxsj  = (xj  - Xj 

dys.j  = (y,  - yj ) 

where 

Xj  ,/j  ,Zj  = coordinates  of  desired  data  point  (t) 
Xj,/j,Zj  = coordinates  of  known  data  point  (j). 

The  desired  elevation  value  is  computed  from 

L Zj/XYZ, 

Z = ia , 

' N 

D I.o/XYZ, 

J=1 


where 

XYZ  j = dxsj  + dysj  + C 
C = constant  (presently  set  at  .01). 


The  interpolation  routine  used  the  same  algorithms  for  storage  of  raw 
contour  data  and  outputting  of  derived  profiles.  Only  the  method  for 
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3.3 


Random  Points 


Unlike  the  contour  or  profile  mode,  the  random  point  mode  com- 
putes coefficients  to  represent  an  entire  area  of  interest.  The  ob- 
jective of  this  link  is  to  take  rather  sparse  discrete  point  data  and 
develop  a terrain  function  which  is  adequate  to  produce  an  orthophoto- 
graph. These  discrete  points  can  be  acquired  from  a photogrammetric 
triangulation,  available  ground  survey  data,  or  points  digitized  from 
available  map  sources.  Coordinate  systems  allowed  as  input  are: 
Universal  Space  Rectangular  (USR),  Geographies,  or  Universal  Transverse 
Mercador  (UTM). 

A multiquadric  function  is  used  to  represent  the  terrain  relief 
over  the  entire  area  for  which  the  orthophotograph  is  to  be  produced 
(see  Section  4.3).  Initially,  the  object  space  area  is  then  grided  in- 
to theoretical  rows  and  columns.  The  number  of  rows  and  columns  to  be 
established  is  designated  by  the  user.  The  points  of  intersection  of 
each  row  and  column  is  defined  as  a nodal  point.  Equation  3.0.1  defines 
the  manner  in  which  these  nodal  point  coordinates  are  applied  in  the 
multiquadric  solution.  The  maximum  number  of  designated  nodal  points 
cannot  exceed  70.  Additional  nodal  points  improve  the  accuracy  of  the 
terrain  representation;  however,  if  an  excessive  number  of  nodal  points, 
relative  to  the  available  relief  data,  are  introduced  the  solution  can 
become  ill-conditioned.  After  the  least  squares  solution  for  the  multi - 
quadric  coefficients  has  been  computed  they  are  used  to  derive  "z" 

coordinates  for  the  new  profiles. 
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3.4  Combinations  of  Data  Types 

Combinations  of  multiple  profile  tapes  (profile,  contour,  and 
random  point  data)  can  be  accommodated  using  the  OOPS  software.  With 
this  option,  profiles  generated  from  several  data  sets,  each  of  which 
covers  only  a portion  of  the  total  area,  can  be  combined  to  generate 
one  complete  profile  set.  This  combination  is  accomplished  in  the 
following  manner.  The  user  inputs  the  total  number  of  data  sets  avail- 
able. Each  data  set  is  then  processed  in  the  normal  manner  and  profiles 
are  output  to  a single  sequential  file  in  the  OOPS  model  system.  After 
all  data  sets  have  been  processed  the  combined  sequential  file  is 
searched  for  profiles  with  common  coordinates.  The  "Z"  coordinates  of 
common  points  are  averaged  and  an  output  tape  formatted  for  the  off-line 
orthophoto  printer  is  output. 

Random  data  points  may  be  included  directly  into  either  the 
profile  or  contour  mode.  These  included  points  are  used  to  strengthen 
either  the  profile  or  contour  solution  in  a specific  area.  The  random 
points  entered  in  this  manner  are  input  in  either  a UTM,  USR  or  geo- 
graphic coordinate  frame. 

The  random  points  are  transformed  from  UTM,  USR,  or  Geographic 
coordinates  to  the  coordinate  system  of  the  data  type  being  processed. 
These  additional  points  are  then  combined  with  the  original  data  set 
and  processed  in  the  normal  manner. 
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MATHEMATICAL  CONCEPTS 


The  mathematics  of  the  OOPS  software  system  can  be  divided 
into  three  general  areas:  (1)  establishment  of  coordinate  systems  and 

coordinate  transformations,  (2)  computation  of  the  photographic  foot- 
print and  digital  data  limits,  and  (3)  derivation  of  the  functional 
coefficients  which  represent  the  terrain  surface. 

4. 1 Coordinate  Systems  and  Transformations 

Orientation  data  can  be  input  in  various  coordinate  systems. 
This  allows  for  complete  flexibility  in  the  choice  of  source  imagery 
to  be  rectified  on  the  OOP.  Initially  the  available  digital  data 
source  must  be  related  to  the  photographic  imagery  from  which  the 
orthophotograph  is  to  be  produced.  In  addition,  the  parameters  re- 
lating the  photographic  imagery  to  an  absolute  coordinate  frame  must 
be  transformed  into  the  model  system  of  the  orthophoto  printer. 

The  photographic  imagery  from  which  the  orthophotograph  is  to 
be  produced  can  be  related  to  any  of  several  absolute  coordinate  frames. 
The  photographic  exposure  station  can  be  defined  in  Universal  Space 
Rectangular  (USR  or  geocentric),  geographies,  absolute  AS-11  model  sys- 
tem, or  an  earth  tangent  LSR  system.  The  photographic  attitude  can  be 
expressed  as  angles  in  the  following  systems:  USR,  AS-11  and  Local 

Vertical  (LV).  If  the  user  desires,  an  orientation  matrix  in  either 
USR  or  LV  may  be  input  to  designate  the  photo  attitude.  The  photo- 
graphic positional  rates  must  be  input  in  a local  vertical  system  or 


-24- 


the  AS-11  system.  The  associated  photographic  attitude  rates  are 
input  as  local  vertical  angular  rates,  a local  vertical  rate  matrix 
or  AS-11  angular  rates.  In  order  to  establish  a unique  orthophoto 
model  system  regardless  of  the  input  data  source,  the  projected  prin- 
ciple point  of  the  photograph  is  established  as  the  origin  of  an 
arbitrary  earth  tangent  local  model  system.  The  "y"  axis  of  this  local 
model  system  is  established  as  being  the  projected  positive  "y"  axis  of 
the  output  photo.  The  "x"  axis  is  then  defined  to  form  a right  handed 
local  system.  Figure  4.1.1  is  an  illustration  of  this  earth  tangent 
model  system. 

The  detailed  mathematics  for  each  specific  coordinate  trans- 
formation is  contained  in  the  subroutine  descriptions  presented  in 
Section  7.  Only  the  general  flow  of  the  required  transformations  will 
be  presented  in  this  section.  Each  function  is  identified  with  the 
section  containing  its  documentation  so  the  reader  can  easily  reference 
the  detailed  mathematics. 

Following  is  a description  of  the  transformations  required 
depending  on  the  type  of  orientation  information  with  the  photographic 
imagery.  Initially  all  orientation  data  is  taken  to  a local  vertical 
system.  Note:  The  notation  f refers  to  a function  for  which  the 
mathematics  is  developed  in  the  subroutine  documentation  (Section  7). 
The  subscript  "n"  refers  to  the  sub-section  that  specifically  documents 
that  function. 
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4.1.1 


USR  Orientation 


<P 

A. 

H 


^7,97  0^  N ' 


a Eo 


T = rt  t 

*LV  ■'US'^Lv  'uSR' 


where 


X 

Y 

Z 

<P 

X 

H 


- are  the  USR  coordinates  of  the  photo  exposure  station 


USR 


a EO 


' US  R 
^U9  R L ' 

Tlv 


- are  the  geographic  coordinates  of  the  exposure  station 


is  the  USR  orientation  matrix 

is  the  rotation  matrix  relating  the  USR  and  local 
vertical  systems 

is  the  local  vertical  orientation  matrix. 


4.1.2  AS-11  Orientation 

The  AS-11  absolute  model  orientation  is  given  as: 

A,  A3  0 
-A3  A,  0 
0 0 1 

Rw  ft 
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^x,  - X " 

Y,  -Y, 

II 

Y,  -Y, 

B B 



LE*  -E,  ^ 

E,  - E„ 

where 


A3 


X.,Y.,E 

X.,Y.,E 


s 

s 

a 

a 


the  model  to  ground  scale  factor 
cos  (AZ) 
sin  (AZ) 

ground  coordinates 

point  of  tangency  of  the  local  model  system 

AS-11  model  coordinates 

center  of  gravity  of  model  system. 


Presently  the  OOPS  software  system  assumes  that  the  local 
model  point  of  tangency  is  given  in  UTM  coordinates.  This  can  be  easily 
expanded  if  additional  user  requirements  are  identified.  This  absolute 
model  orientation  is  converted  to  local  vertical  in  the  following  manner. 


E7.92  (Xg  »Yj 


f, 


7.  S9 


X 

'~BX" 

Y 
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BY 

Z 

1 S R 
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The  azimuth  of  the  local  system  is  defined  as: 

AZ  = arc  tan  (A2,A1). 

The  rotation  to  a local  space  system  is  then  defined  in  functional 

form  as: 

^US  R L S R — ^7.9  1 /^O  / AZ)  . 

The  USR  coordinates  of  the  exposure  station  are  computed  as; 


"x~ 

'^1 

Y 

^US  R L S R 

Y + 

‘ 

USR 

Jlsr 

The  USR  exposure  station  coordinates  are  then  converted  to 


geographies 


The  AS-11  orientation  matrix  is  formed,  rotated  to  the  USR 
coordinate  system  and  then  to  local  vertical 

"^USR  ~ ^USRLSR  ^As-11 
Tlv  ~ ^USRLV  ’"^USR* 

4.1.3  Local  Vertical  Orientation 

If  the  orientation  data  is  entered  in  local  vertical  no  trans- 
formations are  required.  i 

4.1.4  Conversion  of  Orientation  to  Model  System 

After  the  relationship  of  the  photograph  to  the  local  vertical 
coordinate  frame  has  been  established,  the  required  orthophoto  at-ea  is 
projected  to  the  earth's  surface  (see  Section  7.63).  This  defines  the 
footprint  of  the  photograph  on  the  earth's  surface.  This  information 

« 

is  stored  for  later  determination  of  the  orthophoto  limits.  The  photo- 
graphic orientation  as  defined  in  the  local  vertical  coordinate  frame 
is  converted  to  an  AS-11  model  system  in  order  to  output  the  orthophoto 
printer  shutdown  tape.  That  relationship  is  computed  in  the  following 
manner. 

The  projected  principle  point  of  the  photograph  is  established 
as  the  AS-11  model  origin.  The  projected  +Y  photo  axis  is  defined  as 
+Y  model  axis. 


i 

A 
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The  local  vertical  orientation  matrix  is  rotated  to  USR  and 
then  to  a local  space  rectangular  system. 

"^LSR  ~ ^USRLsR  ^USRLV  '"^LV  > 

where 

Tlv  “ Local  vertical  orientation  matrix 
Rlvusr  = Rotation  from  local  vertical  to  USR 
■^LsRusR  = Rotation  from  USR  to  Local  Space  Rectangle. 

The  AS-11  angles  are  extracted  from  the  T^j^  matrix  (see 
Section  7.28). 

The  exposure  station  coordinates  in  the  model  system  are  computed 
as: 


~X^ 

X u S R ^0  A 

Y 

1/^n  R L S R 

Yusr  “ Yqr 

Z 

AS-11 

USR  ^0  _ 

where 

P^YZ]usr  USR  exposure  station 
[XYZIqp  = local  model  origin. 


All  the  AS-11  values  are  output  on  paper  tape  for  use  in  orthophoto 
orientation  of  the  photograph  when  generating  an  orthophoto  on  the  OOP. 


4.2 


Computation  of  Photographic  Footprints  and  Digital  Data  Limits 


One  of  the  required  inputs  to  the  OOPS  system  is  the  coor- 
dinate limits  of  the  digital  data  tape  in  both  the  digital  data  coor- 
dinate system  and  a ground  coordinate  frame  (USR,  GEOGRAPHIC,  UTM). 
Given  this  information,  a translation  and  rotation  representing  the 
relationship  from  ground  to  digital  data  can  be  computed. 

The  ground  coordinates  are  initially  transformed  into  a USR 
coordinate  system  and  then  transformed  to  local  coordinates. 


The  transformation  from  UTM  coordinates  to  USR  coordinates  is 
computed  by: 


^ " ^7.93  (UTM  coordinates) 


^ ^7.5g  (*P  /H)  • 


If  the  ground  coordinates  are  input  in  geographies,  the 
transformation  to  USR  is  simply 


" ^7.59 

— _ USR 


The  transformation  to  LSR  coordinates  is  then  computed  by; 


"x“ 

"x  X, 

Y 

= RT 

‘'USR  Lv  * 

Y - Y, 

Z 

L 3 R 

USR  ^ USR 

Where  the  USR  origin  is  the  projected  principal  point  of  the 
photograph  as  described  in  Section  4.1. 

A transformation  is  then  computed  relating  the  digital  data 
corners  to  the  derived  local  ground  coordinates.  The  equation  used 
for  this  transformation  is; 


— — 

— 

— “ 

— — 

X 

= S 

COSO! 

sin  CL 

X 

- 

X 

Y 

L S R 

-sin  CL 

cos  a 

J 

y 

0 D 

Y 

The  least  squares  solution  for  this  transformation  is  then 
derived  in  the  following  manner; 

A = N“^C 

where 

N = B^B 
and 

C = B^e  . 
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The  B matrix  is  computed  as: 


5f,/6a 

bf,/bo(  bfy/bxbff  ^fy/^Xoff  ^fy/2» 

where 

= (cosa  • XpB  . s+  sin  a . /bb  • s - X(,,f)  - Xlsr 
fy  =(-sin  a • Xpo  • s+  cosa  * Xdd  * * " Xjff  ) " 


The  constant  vector  is  computed  as; 


[E]  = 


Xl  SR  - cosa  • X • s - sin  a • y • s + , 

Yls,  + sina  • X • s - cosa  • y • s+  Xd  f f J 


Using  the  coordinates  of  the  corners  of  the  film  format,  USR 
coordinates  are  projected  ontheearth  (see  Section  7.63).  Knowing 
these  USR  coordinates  and  the  rotation  and  translation  computed  above, 
the  model  format  is  transformed  into  the  digital  data  system  (see 
Section  7.6).  These  model  format  coordinates  in  the  digital  data  system 
are  compared  with  the  limits  of  the  digital  data  tape.  This  comparison 
yields  the  location  of  the  desired  area.  These  computed  boundary  values 
are  used  in  deriving  the  final  output  profiles. 
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4.3 


Multiquadric  Representation 


The  basic  equation  used  in  the  OOPS  multiquadric  mathematics 
is 

Zj  = S + +C)~^ , 


Xj,Yj,Zj  = coordinates  of  the  jth  data  point 

Qj  = coefficient  associated  with  the  ith  nodal 

point 

Xj  ,Yj  = coordinates  of  the  ith  nodal  point 

C = an  arbitrary  constant  term  representing  the 
area  of  influence  of  nodal  points 


The  least  squares  approach  is  used  in  computing  the  value  of 
a,  in  the  above  equation.  The  least  squares  solution  involves  solving 
the  equation 


a = N C 

(n,i)  (n,n)(n,i) 
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N = L B'j  Bj 

(n,n)  1=1  (n,i)(i,n) 

C = £ B^j  ej 

(n,i)  1=1  (n,i)(i,i) 

®1  “ ^®1  1 ^1  1 • * * 3 ^ 


1 

Bi3  = 1.0/(Xj-X,f  + (Y,-Y,)"  +C)^  . 

(1/1) 

The  least  squares  solution  is  applied  differently  to  each  data 
type  due  to  the  different  amount  and  distribution  of  data  and  different 
storage  techniques  used.  The  mathematics  as  explained  above  are  the 
same  for  all  data  types,  the  method  of  formation  is  the  only  difference 
(see  Section  3). 

Profile  generation  is  accomplished  by  evaluating  the  multi- 
quadric equation  for  height  z given  various  x and  y values.  These 
values  are  then  output  to  a magnetic  tape  for  conversion  to  a form 
acceptable  on  the  OOP,  Section  3 explains  profile  generation  for  each 
data  type. 
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5. 


USER’S  INPUT  DESCRIPTION 


The  OOPS  system  requires  that  some  data  be  input  as  card 
images.  These  card  images  can  be  input  either  through  a card  deck  or 
a disc  file.  Contained  on  these  images  is  information  describing  the 
photograph,  input  terrain  data  and  final  model. 

Data  required  to  describe  the  photograph  includes  the  ex- 
posure station  and  orientation.  The  exposure  station  coordinates  can 
be  input  in  geographic,  LSR,  USR  or  AS-11  coordinates.  The  orientation 
coordinate  system  can  be  in  LSR,  USR,  Local  Vertical  or  AS-11.  Other 
information  needed  to  describe  the  photograph  is  photo  type  and  camera 
parameters.  Coordinates  of  the  input  terrain  data  corners  and  the 
associated  ground  coordinates  of  those  points  are  input  to  allow  the 
computation  of  a relationship  between  the  input  terrain  data  and  the 
final  model  system.  Other  information  required  for  terrain  data  pro- 
cessing includes  file  storage  sizes  and  nodal  point  patterns. 

The  model  corners  and  profile  increments  must  be  input  and 
used  for  describing  the  final  model.  This  information  along  with  data 
describing  the  photograph  will  be  used  to  compute  the  information 
needed  to  set-up  the  OOP  model. 

A detailed  description  of  the  setup  of  the  card  images  used 
to  input  all  necessary  data  is  contained  in  the  Off-line  Orthoprinter 
Software  (OOPS)  Computer  Program  Documentation. 
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6. 


SYSTEM  AND  SUBROUTINE  FLOW  DIAGRAMS 


A diagram  illustrating  the  complete  OOPS  system  is  presented 
in  this  section.  This  system  diagram  (Figure  6.1)  illustrates  the 
tree  structure  used  in  overlaying  the  OOPS  software.  Presented  on 
this  diagram  is  the  object  module  name  along  with  the  name  of  the 
individual  subroutines  which  make-up  each  object  module. 

The  tree  structure  is  arranged  such  that  subroutines  which  are 
repeatedly  called  by  different  object  modules  are  contained  in  the  ob- 
ject module  called  MAIN.  During  the  execution  of  the  OOPS  system,  the 
first  object  module  to  be  called  from  the  main  program  is  the  INPUT 
module.  By  exercising  various  options,  the  three  modules  on  the  lower 
portion  of  the  tree,  below  INPUT,  will  be  utilized.  After  completing 
the  computations  in  the  INPUT  module,  program  flow  will  proceed  to 
either  EXTRA,  PROFIL  or  CONTR,  depending  on  the  type  of  data  to  be  pro- 
cessed. Following  the  processing  of  the  data  type,  the  MERGE  module 
will  create  a final  output  tape. 

Subroutine  flow  charts  of  the  major  subroutines  are  presented 
in  the  Off-line  Orthoprinter  Software  (OOPS)  Computer  Program  Documen- 
tation. These  flow  charts  give  overall  logical  flow  and  not  the  de- 
tailed program  structure. 
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7. 


SUBROUTINE  DESCRIPTIONS 


Each  subroutine  developed  for  the  OOPS  software  system  is 
explained  in  this  section.  The  purpose,  entry  points,  calling  ele- 
ments, common  area  used,  subroutines  called  and  link  containing  the 
subroutine  is  all  presented  on  the  following  pages.  A mathematical 
description  is  included  for  subroutines  containing  major  mathematical 
techniques.  Subroutines  involving  storage  or  indexing,  or  those  con- 
sidered standard  contain  no  mathematical  technique  section. 
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7.1  Subroutine  ADD  QAP) 

Purpose:  Determine  location  of  coefficient  in  random  file. 

Entry  Points-  None 

Elements:  I/O  Type  Vorloble  Dimension  Description 

0 I lAP  - Location  on  disc 

Common  Areas:  /PPINT/,/STORE/,/CONST/,/tPDLPR/, 

Subroutines  Used:  None 

Link:  WORK 


7.2  Subroutine  ADDATA 

Purpose:  To  transform  survey  or  terrain  data  from  USR  coordinates 

to  digital  data  coordirKites. 

Entry  Points:  None 

Elements:  None 

Common  Areas:  /UNITS/,/LVTRAN/,/PH0T0/ 

Subroutines  Used:  MATMPY,MOVE,PLHXYZ 

Link:  COR 

Mathematics: 

Given  the  origin  in  Geographic  coordinates,  the  USR  coordinates 

ore  computed  by: 


(7.2.1) 


- — 

<p 

Yo 

“ ^7.60 

X 

3- 

USR 

h 

SCO 


Knowing  the  USR  coordinates  of  the  origin  and  the  USR 
coordinates  of  the  survey  or  terrain  data  points.  The  digital  data  coordinates  ore 
computed  by  computing  the  local  USR  coordinates  from: 


X 

"x*  -Xo' 

Y 

= 

Y*  - Yo 

local 

rJ 

1 

« 

N 

J 

(7.2.2) 


where 


Xj  ,Y,  ,Zj,  = USR  coordinates  of  the  additional  survey  or 
terrain  points. 
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Knowing  the  local  USR  coordinates,  the  digital  data  coordinates 
are  simply  a rotation  and  scale  application: 

(7.2.3) 


where 


*00 
/oo 
^0  0 


cosa  -sin  a 0 
sin  a cosa  0 

0 0 1 


X 

Y 

Z 


o/o 


LOCAL 

USD 


® = ozimoth  ground  to  model  axis 

^eo/a  ~ scale  (digital  data  to  ground) 


Equation  (7.2.2)  ond  (7.2.3)  ore  repeated  for  all  added  points. 


7.3  Subroutine  ADORES  (lAP) 


Purpose: 

Entry  Points: 
Elements: 

Common  Areas: 
Subroutines  Used: 
Link: 


To  determine  location  of  coefficient  in  random  file. 
None 

I/O  Type  Vorioble  Dimension  Description 
0 I lAP  - Location  on  disc 

Apint/,/Store/,/const/,Apdlpr/ 

None 

CREATE 
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7.4  Subroutine  AMXMN  (XMMN,YX,YM) 

Purpose:  To  locate  the  minimum  and  maximum  X and  Y values  of  the 

tope. 

Entry  Points:  None 


Elements: 

I/O 

Type 

Variable 

Dimension 

Description 

1 

R 

XMMN 

(2,4) 

Corner  coordinates 

0 

R 

YX 

(2) 

Maximum  X,  Maximum  Y 

0 

R 

YM 

(2) 

Minimum  X,  Minimum  Y 

Common  Areas:  None 

Subroutines  Used:  None 

Link:  COR 


7.5  Subroutine  ABATES  (XC,ANGLV,ANGR) 

Purpose:  Compute  angular  rates  for  AS  11  shutdown  tape. 


Entry  Points: 

None 

Elements: 

I/O  T>pe 

Variable 

Dimension 

Description 

1 R 

xc 

(7) 

Exposure  stotion  coordinates 

arxi  angles 

1 R 

ANGLV 

(6) 

L.V.  ongle  and  rote 

« 

0 R 

ANGR 

(3) 

AS  1 1 rote  angle 

Common  Areas: 

None 

Subroutirres  Used: 

AS  llAN, CLEAR, GEOPM,MATMPY,MINV,OMATLV,USRLVR,XYZPLH 

Link: 

SDT 
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Mathematics: 


Given  the  Locai  vertical  angles  roll  (a>),  pitch  (fp),  and 
yaw  (a)  and  their  associated  local  vertical  angular  rates  (dw/dt),  (dp/dt),  and 
(da/dt)  this  subroutine  transforms  the  known  vertical  rates  to  the  AS-11  system. 

The  AS-11  orientation  matrix  is  defined  os: 


where 


let 


then 


or 


T*,- 


11 


'Tl»(RusHL*VO^USIIAS-Ii) 


(7.5.1) 


T = local  vertical  orientation  matrix 

IT  “ rotation  matrix  from  USR  to  local  vertical 

Ry,n  AS-ii  ~ '^ototion  matrix  relating  USR  to  the  AS-11 
local  space  rectangular  system 

1^  “ t»)  (^us«**-ii)  (7.5.2) 

^At-ii  ~ I'ur  ^ 

^ ® Vf  * Ta»-ix 
dT^,_ii  = TtydR  + dT^^R 


*^*s-ii  *^L*  s-ii 


(7.5.3) 


where 

A(l,l)  = 


A = <lT^yTjy  from  equation  (7.5.3)  can  be  derived  as: 


0 ' -cos^cos  O)  dot/dt  + sinwd<p/dt  | cos^sin  U)  da/dt  + cos  lo  cRo/dt 
COSTCO*  «da/dt  - sin  w ckp/dt  [ 0 • -sin<p  da/dt  + dw/dt  + dw/dt 

I — — — — 

-ein  OJ  cos<p  da/dt  - cosoj  ckp/dti  sin  <p  da/dt  - dco/dt  i 0 


(7.5.4) 


dRc(Ryy^l^y)  d Ry  j , ^ j j + dRy 


JSA  L« 


Ryj 


A AS-11 
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The  derivoHve  of  rofoMon  matrix  from  USR  fo  rfie  local  space 
lystem  of  the  AS-11  plotter  is  zero  and  this  equation  becomes: 


Since 


*1-11  ~ 0 


dR  = dR, 


us*  L»  * Kgs*  AS-11 


this  equation  can  be  rewritten  os 

dR  = [dT^  • ♦ dT^  • i -p 

^ V V V V V 


(7.5.5) 

(7.5.6) 


where 


cosy''  -sin  0 

sin  y"  cosy  0 

0 0 1 


J 


1 0 0 
0 sinrp*^  cos?)' 

0 -cos?)'  sin?)' 

-sinX'  cosX'  0 
-cosX'  -sin  X'  0 
0 0 1 


dT’  , dT^  , and  dT^  are  derivatives  of  these  rrKJtrices  with 

y X' 

respect  to  time. 

The  AS-11  rates  can  then  be  extracted  from  the  rate  matrix 
derived  using  equatiori  (7.5.3). 
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d^/dtAS-11  =dT*,_„(1,3)/cos<p 

doj/dt  AS-11  = (T(2,3)  + sin  w sin  <pd<p)/(cos  w coscp) 

dK/df  AS-11  = (T(l,2)  + sin  <p  cos  X d<p)/(-cos^sin  x) 


where 

tp,tj3,x  are  the  AS-11  orientation  angles 
TO,  3)  = cos^d^ 

T(2,3)  = cosojcoscpdoj- sinO)sin(pd<p 
TO, 2)  = -sin^cosxd((5  - coscpsin  xd  X 

These  derived  rates  are  a function  of  time  and  must  be  converted 
to  bed  function  of  the  scan  angle  for  use  in  the  AS-11  system.  Thus: 

d^/da  = &p/d\  • dt/da 
6w/iHol  = dw/dt  • dt/da 
dx/da  = dx/dt  • dt/da. 

7.6  Subroutine  ARXY  (XY) 

Purpose:  To  find  maximum  and  minimum  corner  values  needed  for 

stripping  tape. 

Entry  Points:  None 

Elements:  I/O  Type  Vorioble  Dimension  Description 

I R XY  (2,4)  Image  comer  measurement 

Common  Areas:  APDLPR/,/UNITS/,/GR0UND/,/0PTI0N/,/PH0T0/, 

/MINMAX/,/MODELP/,/LVTRAN/ 

Subroutines  Used:  AMXMN  ,MATMPY,PLHXYZ  ,MOVE 

Link:  COR 
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Mathematics: 


Given  the  origin  in  geographic  coordinates,  the  USR 
coordinates  are  cornputed  by: 


Yo 

“ ^7.60 

(7.6.1) 


Knowing  the  USR  coordinates,  the  rotation  matrix  relating  USR  to  LV  is 


“ ^7, 


83 


■o 

Zo 


(7.6.2) 


Given  the  projected  film  comers  in  geographic  coordinates, 
the  USR  coordinates  are  computed  using  Equation  (7.6.1).  The  Local  Vertical 
coordinates  are  then  computed  by: 


Xt  y 

Yt, 

Y,  -Yo 

where 

X,  ,Yj  ,Zj  = USR  coordinates  of  comer  i 

= Local  Vertical  coordinates  of  comer  i 


The  digital  data  coordinates  of  the  projected  film  comers  are  then  computed  by: 


where 


\>o 

Xuv  ” Xof r 

ybo 

«=  R 

Yu  * - Yq  p f 

^0 

„Zuy  - Zcrr_ 

(7.6.4) 


R «»  rotation  matrix  (Section  7.88) 

^ translation  elements  (Section  7.88) 
A " scale  relating  model  to  ground. 
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Comparisons  ore  then  made  between  the  digital  data  coordinates 
of  the  projected  film  format  and  the  input  limits  of  the  digital  data  tape.  The 
minimum  ond  maximum  values  of  these  comparisons  are  used  for  extroction  of  data 
and  final  computotions.  The  area  covered  by  this  final  output  tape,  in  model  scale, 
con  be  computed  by  transforming  the  minimum  arxJ  maximum  values  from  digital  data 
to  ground  and  then  ground  to  model.  These  transformations  ore  performed  as  follows: 


where 


- |Rm/« 


.X  - 


Yorr 


(7.6.5) 


where 


- ...  » • r cos(AZ)  sin  (AZ) 

[-sin(AZ)  cos(AZ) 

AZ  = azimuth  relation  model  to  ground 

x,y  = digital  data  coordinates 

Xg,Yg  = digital  data  corner  in  local  vertical 
ground  coordinate  system 

X = scale  digital  data  to  ground 

Yoff  = ground  translation  volues. 


cos(AZIM)  -sin  (AZIM) 
sinC^ZIM)  cos  (AZIM) 


(7.6.6) 


AZIM  = azimuth  relating  model  to  ground 
= scale  relating  ground  to  model 
X,  ,y,  = model  coordinate  of  digital  dato  corners. 


7.7  Subroutine  ASTRP 


Purpose: 

Entry  Points: 
Elements: 

Common  Areas: 
Subroutines  Used: 
Link: 


To  compute  the  delta  X between  profiles. 

None 

None 

Apdlpr/,Ainits/,/const/ 

None 

STORE 


-48- 


7.8 


Purpose: 

To  compute  the  ASH  angles  from  an  LSR  matrix. 

Entry  Points: 

None 

Elements: 

I/O  Type  Variable  Dimension 

Description 

1 R TLSR  (3,3) 

LSR  matrix 

0 R ANASn  (3) 

ASH  angles 

Common  Areas; 

None 

Subroutines  Used: 

None 

Link: 

SDT 

Subroutine  AS  11M 

(ANG,OMAT) 

Purpose: 

To  compute  photo  to  ground  matrix 

for  ASH. 

Entry  Points: 

None 

Elements: 

I/O  Type  Variable  Dimension 

Description 

1 R ANG  (3) 

ASH  angles  (x,( 

0 R OMAT  (3,3) 

AS  11  matrix 

Common  Areas: 

None 

Subroutines  Used: 

None 

Link: 

Mathematics: 

INPUT 

Given  the  AS-11  angles  the  matrix  relating  model  coordinates 
to  i^iotograph  coordinates  is  defined  as; 

(7.9.1) 

cos X cos^  - sin  x sin  (O  sin  <p  sin  X cos^>  + cos X sin  u sin  <p  -cos u sin  <p 
" »ln  X cos  0)  eos  X cos  u sin  U) 

^cosx  sin  ^ -t-sin  X sin  («j  cosfp  sin  X sin  ^ - cos  x sin  (Ocos^  coswcos^ 


7.10  Block  Data 


Purpose: 

Entry  Points: 

Eletnents: 

Common  Areas: 

Subroutines  Used: 

Link: 

7.11  Subroutine  BMAT  (X,Y,B) 

Purpose:  Com(.Hjte  the  B matrix. 

Entry  Points:  None 

Elements:  I/O  Type  Vgrioble  Dimension  Description 

I R X - X value 

I R Y - Y value 

OR  B (70)  B matrix 

Common  Areas:  /I'JODES/j/CONST/ 

Subroutines  Used:  None 

Link:  WORK 

Mathematics: 

Computes  the  matrix  of  the  partial  derivatives  over  all  nodes. 
This  matrix  is  defined  as: 

i 

BO)  = 1.0/CAX3  + AY3 +C)  (7.11.1) 

where 

AX  = X - XNODE(I) 

AY  = y - YNODE(I) 

C = constant  value  related  to  data  amount 
X,Y  = digital  data  coordinates 
XNODE(l),YNODE(l)  = coordinates  of  node  I . 

I = 1,  number  of  nodes 


To  set  up  common  blocks  with  associated  fixed  no. 

None 

None 

Ajnits/,/idxbi/,/index/,/modelp/ 

None 

MAIN 
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7.12  SubrooHne  BMATRX  (X,Y,B) 

Purpose;  To  compufe  the  B matrix. 

Erttry  Points:  None 

Elements:  I/O  Type  Vorioble  Dimension  Description 

I R X -•  X value 

I R Y - Y value 

OR  B (70)  B matrix 

Common  Areas:  /SlODES/j/CONST/ 

Subroutined  Used:  None 

Link:  CREATE 

Mothemotics:  See  BMAT  (Section  7.11). 


7.13  Subroutine  BMATX  (X,Y,B) 

Purpose: 

Entry  Points: 

Elements: 


Common  Areas:  /1slODES/,/CONST/ 

Subroutines  Used:  None 

Link:  EXTRA 

Mathematics:  See  BMAT  (Section  7.11). 


Form  the  B matrix. 

None 

I/O  Type  Variable  Dimension  Description 

I R X - X model  value 

I R Y - Y model  value 

OR  B (70)  B motrix  computed  by  the 

contribution  of  input  X&Y 
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7,14  Subroutine  BOUND 


Purpose: 

Entry  Points: 
Elements: 

Common  Areas: 
Subroutines  Used: 
Link: 


To  read  input  tope,  select  necessary  dato  and  store  to 
random  disk  file. 

None 

None 

/Store/,/Tpdlpr/,/data/,/ppint/,Ajnits/,/const/ 

BUFFIN,DISC,MXMN,RWRITE,XFILL,YFILL,MOVE,PK 

STORE 


7.15  Subroutine  BUFFIN  (IFLAG,NWD) 


Purpose: 
Entry  Points: 
Elements: 


Common  Areas: 
Subroutines  Used: 
Link: 


To  read  in  data  from  tapes  created  cn  18-bit  machines. 

None 

I/O  Type  Vorioble  Dimension  Description 

I I IFLAG  - Flag  for  contour  or  profile 

I I NWD  - Number  of  16-bit  words 

A>ata/ 

MTREAD 

STORE,CONTR 


7.16  Subroutine  BUFFOT 


Purpose: 

Entry  Points: 
Elements: 

Common  Areas: 
Subroutines  Used: 
Link: 


To  bring  in  2 16-bit  words  and  create  1 18-bit  word. 

None 

None 

/DATA/,/BUFFER/ 

MTREAD 

MERGE 
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7.17  Subroutine  CLEAR  (A,N) 


Purpose: 
Entry  Points: 
Elements: 


To  set  N elements  of  array  A to  zero 
None 


I/O  Type  Variable 

I/O  R A 
I I N 


Dimension  Description 

(1)  Array  to  be  set  to  zero 

Number  of  words  in  A to  be 
set  to  zero 


Common  Areas;  None 

Subroutines  Used:  None 

Link:  MAIN 


7.18  Subroutine  CONTF  (Z) 


Purpose: 

Entry  Points: 
Elements: 

Common  Areas: 
Subroutines  Used: 
Link; 

Mathematics; 


Will  compute  a Z value  given  on  X and  Y coordinate. 

None 

I/O  Type  Voriable  Dimension  Description 
0 R Z - Computed  Z value 

/STORE/,/CONST/,/t>ATA/,/PPDLPR/,/UNITS/,/PARMS/ 
ADDRES,BMATRX,UNPRPK 
CREATE 


Given  the  x.y  coordinates,  the  storage  address  for  the 
coefficients  con  be  computed  (Section  7,3).  The  center  of  this  coefficient  area  is 
subtracted  from  the  coordinates  to  put  them  in  a local  system. 


where 


X,  (7.18.1) 

Xt=Y-Y„,  (7.18.2) 

X,Y  = input  digital  data  coordinates 
,Yqp,  = coordinates  of  center  of  coefficient  area. 
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Form  B mofrix  where 


[B]  = f,.„  (X,,YJ. 
Compute  the  Z coordinate  from: 

Z = [B]  • [P] 

where 


(Section  7.49). 


(7.18.3) 

(7.18.4) 


7.19  Subroutine  CONTR 


Purpose: 

Entry  Points: 
Elements: 
Common  Areos; 

Subroutines  Used: 
Link: 


To  locate  contour  header  from  digital  data  tape  and 
begin  collecting  data. 

None 

None 

/index/,/modelp/,/addres/,Ajnits/,/data/,Ahoto/, 

/PARAM/ 

IMVSTR,IRANIO,ISTOR,LSTOR,MVSTR,SORTBL,BUFFIN, SORTER 
CONTR 


7 JO  Subroutine  CORNER 


Purpose: 

Entry  Points: 
Elements: 

Common  Areas: 
Subroutines  Used: 
Link: 


Determine  minimum  and  maximum  values  of  model  system. 

None 

None 

APD  LPR/,/MI  NMAX/,/PH  OTO/ 

MOVE 

EXTRA 
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7.21 


Subroutine  CORUSR  (ICARD,ITYPE) 


Purpose; 


£ntry  Points: 
Elements: 


Common  Areas; 
Subroutines  Used; 
Link: 


Reads  digital  data  corner  information  in  ii  age  coordinates 
end  object  coordinates.  A matrix  is  then  formed  to  relate 
these  two  coordinate  systems. 

None 

I/O  Type  Varioble  Dimension  Description 

I I ICARD  - Option:  data  input  (CARDS, TAPE) 

I I ITYPE  - Option:  data  input  (USR,GEC,UTM) 

/UNITS/,/t)lGlT/,/OPTION/ 

DMSRAD,PLHXY2,TRANS,ADDATA,ARXY 

INPUT 


Mathematics: 


Given  the  UTM  coordinates  of  a point,  the  USR  coordinates 
•ore  computed  by  first  converting  to  geographic  coordinates.  This  computation  is 
computed  by; 


where 


(7.21.1) 


XyT*(/Y„TH  = UTM  coordinates  of  the  corner. 


is: 


KrMwing  the  geographic  coordinates,  the  conversion  to  USR 


(7.21.2) 
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7.22  Sobfoutine  DDTOMD(XT,YT,XPNT,YPNT) 


Purpose: 
Entry  Points: 
Elements: 


Common  Areas: 
Subroutines  Used: 
Link; 


To  convert  digitol  data  to  model  coordinates. 
None 


I/O 


0 

0 


Type  Voriable 

R XT 

R YT 


Dimension 


R 

R 


XPNT 

YPNT 


Avtran/,Ahoto/,/param/ 

None 

CREATE 


Description 

Input  digital  data 
coordinotes 


Model  coordinates 


Matbemutics: 

Given  the  digital  data  coordinates  of  a point,  the  model 
coordinates  are  computed  by  converting  the  digital  data  to  Universal  Space 
Rectangular  (USR)  coordinates  and  then  tronsforming  these  USR  coordinates  to 
model  coordinates. 


by: 


where 


The  tronsformation  of  digital  data  to  USR  is  accomplished 


“■  — 

- 

““ 

- — 

r-  — 

X 

= 

cos  Of 

sin  a 

X 

• ^s/oB 

Xb 

Y 

-sin  a 

cos 

y_ 

-< 

o 

\ 

a = azimuth  relating  digital  data  to  ground 
x,y  = digital  data  coordinates 
Xp,Yj  = origin  of  digital  data  in  USR  coordinates 
^»/oo  ^ scale  (ground  to  digital  data). 
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The  conversion  from  ground  (USR)  to  model  is  computed  by 
the  applying  of  scale  and  a rotation  matrix: 


where 


cosy  -sin  y 
sin  y cosy 


X.,Y,  = model  coordinates 

y = azimuth  relating  ground  to  model 

X = scale  relating  ground  to  model 

X,Y  = USR  coordinates. 


(7.22.2) 


7.23  Function  DISC  (L,X,YLM,YLIM) 


Purpose: 
Entry  Points: 
Elements: 


To  compute  location  stored  in  random  file 
None 

I/O  Type  Vorioble  Dimension 

III 


Common  Areas: 
Subroutines  Used: 
Link: 

Mathematics: 


I R X 

I R XLM 

I R YUM 

APDLPR/,/EPINT/,AJNITS/ 

None 

PRCFIL 


Description 

Profile  number  within  1 column 
X value  of  profile 
Upper  limit  of  block 
Lower  limit  of  block 


Compute  the  disc  storage  location  where 

DISC  = IBASE  + N 

IBASE=  10*((IX  - IXMIN)/IBLOCK)*((IYMAX  - IYMIN)/lBLOCK) 
N = (L-  1)  * ((lYMAX  - IYMIN)/IBLOCK)  + IPROF  + 1 


(7.23.1) 

(7.23.2) 

(7.23.3) 
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r 


where 


IX  = X value  of  profile 
IXMIN  = minimum  x value 
IBLOCK  = block  size 
lYMAX  = maximum  y value 
lYMIN  = minimum  y value 
L = profile  number 

IPROF  = blocks  previously  stored  from  profile  L. 


7.24  Subroutine  DMSRAD  (L,M,R,B) 


Purpose: 
Entry  Points: 
Elements: 


To  convert  deg,,  min.,  and  sec,  to  radians. 


I/O  Type  Variable  Dimension  Description 

III  ~ Integer  degrees 

I I A/I  - Integer  minutes 

I R R Seconds  (real) 

Orb  - Radians 


Common  Areas:  None 

Subroutines  Used:  None 

Link;  MAIN 


7.25  Subroutine  DOMLV  (ANG.DOM) 


Purpose: 
Entry  Points: 
Elements: 


To  compute  the  local  vertical  rate  matrix 
None 

I/O  Type  Voriable  Dimension  Description 

I R ANG  (6)  L.V . ongles  ond  rates 


Common  Areas:  None 

Subroutines  Used:  None 

Link:  SDT 


L.V.  rate  matrix 


J 


Mathematics: 


(7.25.1) 


A(1J) 

A(l,2) 

A(l,3) 

A(2,l) 

A(2,2) 

A(2,3) 

A(3,l) 

A(3,2) 

A(3,3) 


cos<p  cosa  da  - sin  ot  sin  dip 
-cos^sin  a da  - cosa  sin  <p  dtfi 
cosifi  dip 

-COSO)  dw  sin  <p  sin  a - sin  OJ  cos<p  cosa  dp  + sin  W sin  p sin  a da 
+ sin  o)  cosa  dco  + cosa)  sin  a da 

-cosu  sin  p cosa  doi  - sin  u cosp  cosa  dp  - sin  U)  sin  p sin  a da 

- sin  o)  sin  a du  -t-  cos  a)  cos  a da 

COSO)  cosp  db)  - sin  h)  sin  p dp 

sin  b)  sin  p sin  a db)  - cosb)  cosp  sin  a dp  - cosb)  sin  p cosa  da 
-t-  cosb)  cosa  dbt  - sin  b)  sin  a da 

sin  b)  sin  p cosa'  db)  - cosb)  cosp  cosa  dp  + cosb)  sin  p sin  a cia 

- cos  b)  sin  a db)  - sin  b)  cos  a da 

-sin  b)  cosp  db)  - cosb)  sin  p dp 


7.26  Subroutine  EFIT 


Purpose: 

Entry  Points: 
Elements: 

Common  Areas: 
Subroutines  Used: 
Link: 


Evaluate  the  fit  of  the  coefficients  using  known  coordinates. 

None 

None 

/UNITS/,/PARMS/,/CONST/,/NODES/,/tPINT/ 

BMATX,MATMPY 

EXTRA 
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Mathematics: 

Using  the  x,y  image  coordinates,  the  matrix  of  partial 
derivatives  can  be  computed  from 

«=  f7.i3  0c,y)-  (7.26.1) 

Having  computed  the  B matrix,  the  associated  Z value 
can  be  computed  from 

Zc  = [B]*[P]  (7.26.2) 

where 

[P]  = coefficients  derived  from  normal  solution  of 
multiquadratic. 

Given  the  "Z"  value  read  from  the  input  digital  data  tape 
and  the  newly  computed  "Z"  value,  a difference  can  be  computed 

^Z  = Z|,  - Zj  (7.26.3) 

where 

Z„  *=  "Z"  read  from  tape 
Zj  = computed  Z. 

The  sum  of  the  squares  of  AZ  is  used  to  compute  the 

mean  error  in  Z . 


where 


N = number  of  points  evaluated. 


(7.26.4) 


I 


7.27  Subroutine  EVAL  (7ME) 


Purpose: 

Entry  Points: 
Elements: 

Common  Areas: 

Subroutines  Used: 
Link: 

Mathematics: 


To  evaluate  fit  of  derived  coefficients. 

None 

I/O  Type  Vorioble  Dimension  Description 
0 R 21ME  - Mean  error  in  Z 

/const/,/data/,/^tore/,apint/,AJNits/,/Tpdlpr/, 

/PARMS/ 

RREAD,VAL 

WORK 


Using  the  x,y  image  coordinates,  the  matrix  of  partial 
derivatives  can  be  computed  from: 


= fv.iiOt,/)* 


(7.27.1) 


Having  computed  the  [B]  matrix,  the  associated  Z 
value  can  be  computed  from 

Zj  « [B][P],  (7.27.2) 

where 

[P]  = coefficients  derived  from  normal  solution  of 
multiquadratic. 

Given  the  *Z"  value  read  from  the  input  digital  data  tape 
and  the  newly  computed  "Z"  value,  a difference  can  be  computed. 

AZ  *=  Z,  - Zj  (7.27.3) 

where 

Z,  » "Z"  reod  from  tape 
Zj  » computed  "Z". 
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error  in  Z . 


where 


The  sum  of  the  squares  of  AZ  is  used  to  compute  the  mean 

Z„  = E (AZ.)^AJ  (7.27.4) 

1=1 

N = number  of  points  evoluated. 


7.28  Subroutine  EXLVAN  (ANGLS) 

Purpose:  To  extract  local  vertical  angles  and  rates 

Entry  Points:  None 

Elements:  I/O  Type  Vorioble  Dimension  Description 

0 R ANGLS  (6)  L.V . angles  and  rates 

Common  Areas:  Blank 

Subroutines  Used:  None 

Link:  SDT 

Mathematics; 

Given  the  local  vertical  orientation  matrix  and  local  vertical 
rate  matrix,  the  local  vertical  angles  and  rates  can  be  extracted. 


- ore  fan 


CO  = arc  tan 


O = ore  tan 


rOMQ,3)  T 
[_OM(3,3)/cos  wj 

r OM(2,3)l 

LoM(3,3) j 

rOM(i,1)1 

LomO,'2)J 


(7.28.1) 

(7.28.2) 

(7.28.3) 
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i 


J 


dp  = DOM  (l,3)/co*v> 

dw  “ [DOM  (2,3)  + sin  CO  sin  ^ dp]/(cos  faj  cos(p) 
da  ■ [DOM  (1,2)  + cos  a sin  <pdip]/(-cos(psina) 

where 

r cos^sina  cos^  cosa  sin  <p 


I 


(7.28.4) 

(7.28.5) 

(7.28.6) 


OM 


•eln  CO  sin  <p  sin  a 
-cos  CO  cosa 


•sin  CO  sin  <p  cosa 
+COSCO  sin  a 


sin  CO  cosip 


(7.28.7) 


-cos  CO  sin  psin  a 
■fsin  cocos  a 


-cos  CO  sin  <p  cosa 
-sin  CO  sin  a 


cos  CO  cos 


DOM(l,2)  = -sinipcosackp-sinacosipda  (7.28.8) 

DOM(1,3)  = co$<pd<() 

DOM(2,3)  = costocospdco-sin  vsin  codp. 


7.29  Subroutine  FPBOND 


Pwpote; 

Entry  Points: 
Elements: 
Common  Areas: 
Subroutines  Used: 
Link: 


Arranges  footprint  corners  so  first  point  has  minimum  X coordinates. 
Also  limits  of  boundaries  are  determined. 

None 

None 

Zl=PDLPR/,/MINMAX/ 

MOVE 

STORE 


1 
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7.30  Subroutine  GENODE 


Purpose: 

Entry  Points: 
Elements: 

Common  Areas: 
Subroutines  Used: 
Link: 

Mothemotics: 


To  generate  Hie  coordinates  of  the  nodes. 

None 

None 

/I=pdlpr/,Aiodes/,/const/,Ainits/ 

None 

STORE 

The  coordinates  of  the  nodes  are  computed  from  Hie  equations: 


NODXO)  = DELX*(J-1)  - XSIZE  (7.30.1) 

NODYO)  = DELY*(K-1)  - YSIZE 

where 

NODX(i),NODY(l)  = coordinates  of  node  I 


XSIZE 

= block  size  in  X 

YSIZE 

= block  size  in  Y 

DELX 

= (2.*xsize)Ajodex 

NODEX 

= total  number  of  nodes  in  X 

DELY 

= (2.*YSlZE)/NODEY 

NODEY 

= total  number  of  nodes  in  Y 

J 

= variable  from  1 NODEX 

K 

«=  variable  from  I -*  NODEY 
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7.31  Subroutine  GEOPM  (XYZ,PLH,XYZPLH) 

Purpose:  To  provide  tbe  matrix  of  partial  derivatives  of  geocentric 

X,Y,Z  with  respect  to  geographic  (p,X,H  . 

Entry  Points:  None 


ElemenK; 

\/0 

Type 

Variable 

Dimension 

Description 

1 

R 

XYL 

(3) 

Geocentric  coordinates  X,Y,Z 

1 

R 

PLH 

(3) 

Geographic  coordinates  o,X,H 

0 

R 

XYZPIH 

(3,3) 

Matrix  of  partial  derivatives  of 

X,Y,Z  with  respect  to<p,A,H 

Common  Areas:  /GEODET/ 

Subroutines  Used:  None 

Link:  SDT 

Mathematics: 

Given  the  geocentric  coordinates  and  geographic  coordinates, 
o matrix  of  partial  derivatives  of  geocentric  X,Y,Z  with  respect  to  geographic 
^,X,h.  The  matrix  is  composed  of  the  following  elements: 

(1.1)  = -[A/(l  - e®  sin®  ^ cosX  + [[«?  (X/(l  - e®  sin®  sin<p  coscp] 

/(l-e*sin®<p)]  cos  ((>  cos  X 

(1.2)  - -Y 

(1.3)  = cosipcosX 

(2.1)  = -[(X/(l  - e®  sin®  ^)^)  + h]  sin^sin  X + [[e®  (A/(l -e®  sin®  VJ)^)  sin  Ip  cos<p] 

/(I  -e^  sin®  <p)^  ] coi(p  sin  X 

(2.2)  = X 


(2,3)  “ cos^sinX 

(3,1)  = cos^[(A/(l  - e®  sin®  :p)^)(l -e®)+ h]  + [[e®  [A/(l  - e®  sin®  V5)^]  sin  ip  cos  ip] 


I 


(3.2)  = 0 

(3.3)  B sinv> 
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(7.31.1) 


where 


A •=  semJ-fnojor  axis  of  the  earth 

e?  = eccentricity  squared. 


Subroutine  GEUTM(OPH!,OlAM,PHI1,ALAM,OK,X,Y) 

Purpose: 

To  transform  geographic 

coordinates  to  UTM  coordinates. 

Entry  Points: 

None 

Elements: 

k9 

IlE£ 

Variable 

Dimension 

Description 

1 

R 

OPHl 

- 

UTM  origin  in  O 

1 

R 

OLAM 

- 

UTM  origin  in  X 

1 

R 

PHIl 

- 

Input  <p 

1 

R 

A LAM 

- 

Input  X 

I 

R 

OK 

- 

SCALE 

0 

R 

X 

- 

Computed  UTM  value 

0 

R 

Y 

- 

Computed  UTM  value 

Common  Areas: 

/GEODET/ 

Subroutines  Used:  None 

Link:  CORD 

Mathematics: 


Given  the  geographic  coordinates  of  a point,  the  UTM 
coordinates  can  be  computed  from 

X ='’  OK  (XPR)+  BX 
Y = OK  CYPR) 

where 

OK  = scale 

BX  ^ air  base  component  in  X 

XPR  = EN  + (DLC  + DLC=(l-(tan<Pi)3  + ET)/6 

+ DLC?((5-  18(ton<Pi)®  + (ton<Pj)«)/120) 

EN  = A/(l-ea(sin(Pif)^ 

DLC  = AX  cos<p^ 

ET  •=  Ce»/(l-e=)](cosV>,y' 
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YPR  = S + EN((DLC  - sin  <p^  AX)/2  + DLC=*  sin  AX 
(5  - (fon<pj)3/24) 

S = A(l-e®)[Aj(<Pi -<Po)  - (B/2)(sin(2cPj)  - sin(2(Po)) 
+ C/4  (sin  (4<pi)  - sin  (4<Po  ))  - (D/6)(sin  (&p^) 
-$ln(6p(,))] 

AX  = \,  -Xj 

A = semi-major  axis 

e*  = eccentricity 

(Pq  ,Xo  = geographic  coordinates  of  origin 
^j,Xj  = geographic  coordinates  of  point 

A1  = 1 + 3/4  e=  + 45/64  e*  + 175/256  e« 

B «=  1/2C3/4  «?  + 15/16  e*  + 525/512  e®  ] 

C = 1/4C15/64  e*  + 105/256  e^"  ] 

D = 1/6  [35/512  e®] 


7.33  Subroutine  GNODES  (NODX,NODY) 


Purpose: 
Entry  Points: 
Elements: 


Common  Areas: 
Subroutine  Used: 
Link: 


Compute  nodal  points  for  multi-quad  solution. 

None 

I/O  Type  Vorioble  Dimension  Description 

I I NODX  - Number  of  nodes  in  X 

I I NODY  - Number  of  nodes  in  Y 

Apdlpr/,/Nodes/,/const/,/Units/,/Tpint/ 

None 

EXTRA 
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The  equation  for  the  computation  of  the  nodel  point 


7.34 


coordinates  are: 

NODXO)  = XMIN  + DELX  * (J-1) 

NODY(l)  = YMIN  + DELY  * (K-1) 

where 

NODXO)/NODY(1)  = coordinates  of  node  I 
XMIN  = minimum  x value  for  computation  area 
YMIN  = minimum  y value  for  computation  area 
DELX  = (XMAX  - XMlN)/(NODEX-l) 

XMAX  = maximum  x value 
NODEX=  number  of  nodes  in  x 
DELY  = (YMAX  - YMIN)/(NODEY  - 1) 
YMAX  = maximum  y value 
NODEY=  number  of  nodes  in  Y 
J = variable  from  1 -*NODEX 
K = variable  from  1 -*NODEY 


(7.33.1) 


Subroutine  IMOVE  (1A,1B,N) 

Purpose: 

Integer  move 

of  data  in  labeled  common  block  /DATA/ 

Entry  Points: 

None 

E lements; 

I/O  Type 

Variable  Dimension 

Description 

1 1 

lA 

First  address  of  array  to  move 

1 1 

IB 

First  address  of  array  moving  to 

1 1 

N 

Number  of  words  to  be  relocated 

Common  Areas: 

/t)ATA/ 

Subroutines  Used: 

None 

Link:  MERGE,CONTR 
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7.35  Subroutine  IMVSTR 


Purpose; 

Entry  Points: 
Elements: 
Common  Areas: 

Subroutines  Used: 
Link: 

Motbematics: 


To  transform  boundary  limits  to  new  coordinote  system. 

None 

None 

/modelp/,/index/,Ainits/,/1>aram/,/data/,/lvtran/, 

Apdlpr/,/1’hoto/ 

None 

SORTPR 


Given  the  digital  data  corners  in  the  digital  data  coordinate 
system,  the  conversion  to  USR  coordinates  is  computed  as; 


cos  (AZ) 

sin  (AZ) 

-sin  (AZ) 

cos  (AZ) 

• ^s/dd  " 

J0R8. 

(7.35. 1) 


where 


AZ 

^/oo 


azimuth  reloting  ground  to  digital  data 
image  coordinates  of  point  i 
scale  relating  ground  to  digital  data 
USR  coordinates  of  digital  data  origin 


With  these  newly  computed  USR  coordinates,  the  model 
coordinates  are  derived  from: 


where 


cos  (AZIM) 

-sin  (AZIM) 

Xu 

1 

sin  (AZIM ) 

cos  (AZIM) 

i.oooA,/^, 


AZIM  = azimuth  relating  ground  to  model 
Xj  ,Yj  = USR  coordinates  of  point  i 

= scale  relating  ground  to  model. 


(7.35.2) 
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7.36 


7.37 


A comporison  between  the  four  digital  data  corners  and  the 


coordinates  of  the  model  comers  will  define  the  area  over  which  profiles  are  to  be 
created  and  the  number  of  profiles  to  be  generated. 

Subroutine  INCLUD 

(1C) 

Purpose: 

To  Include  additional  data  points. 

Entry  Points: 

None 

Elements: 

I/O  Type  Variable  Dimension  Description 

1 1 1C  - Number  of  profiles  currently 

being  processed 

Common  Areas: 

/DATA/,/MODELP/,/LASTPT/,/lJNlTS/ 

Subroutines  Used; 

None 

Link; 

DERIVE 

Subroutine  IN  PHOT 

Purpose; 

To  input  data  necessory  for  computing  model  and  photo 
varlobles. 

Entry  Points; 

None 

Elements: 

None 

Common  Areas: 

Ajnits/,/param/,/geodet/,/addpt/,/off/,/option/, 

/Unpar/,/ground/,/photo/ 

Subroutines  Used: 

clear, DMSRAD,EXLVAN>\ATMPY, MOVE,  PUiXYZ,PROJ, 

ptape,rgeo,rusr,usrlvr,xyzplh,omatlv,domlv, 

RAS11,RUTM 

Link: 

INPUT 

Mothematics: 

If  the  exposure  station  coordinates  are  input  in  USR  or  LSR 

coordlrxstes,  o tronsformation  from  USR  to<p,X,h  must  be  performed 

<P 

X 

X 

“ ^7.99 

Y 

h 

Z 
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7.38  Subroutine  IRANIO  (ICUR,ITYPE,IPRES) 


Purpose: 
Entry  Points: 
Elements: 


Common  Areas: 
Subroutines  Used: 
Link: 


Initialize  random  storage  routine. 

None 

I/O  Type  Variable  Dimension  Description 

I I ICUR  - Initialize  storage  areo 

I I ITYPE  - Not  used 

I I IPRES  - Precision  of  word 

AJ  N ITS/,/DA  '•A/,/RAND  B U/,/PDP  I N D/ 

None 

PROFlL,EXTRA,CONTR 


7,39  Subroutine  ISTOR 


Purpose: 

Entry  Points: 
Elements: 

Common  Areas: 
Subroutines  Used: 
Link: 


Set  up  common  areas  and  clear  arrays  to  be  used  in  STOR 
routine. 

None 

None 

/M0DELP/,AJNITS/,/IDXBI/,/ADDRES/,/DATA/ 

KLEAR 

SORTPR 


Subroutine  KLEAR 

(1DATA,N) 

Purpose: 

Clear  integer  array. 

Entry  Points: 

None 

Elements; 

I/O  Type 

Variable 

Dimension 

Description 

I/O  1 

IDATA 

(N) 

Array  to  clear 

1 1 

N 

- 

Number  of  members  of  array  to  cleor 

Common  Areas: 

None 

Subroutines  Used: 

None 

Link: 

CONTR 
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7.41  Subroutine  LIMIT  (X ,YINTR,RMXY) 


Purpose; 


To  compute  maximum  and  minimum  Y value  given  on  X value. 
Also  gives  X and  Y limits  of  the  sheet. 


Entry  Points: 

Nor« 

Elements; 

I/O  Type  Variable 

Dimension 

Description 

1 R X 

- 

X value 

0 R YINTR 

(2) 

Y intersections  given  Y 

Common  Areas; 

0 R RMXY 

/Photo/,/1fpdlpr/ 

(4) 

Max  & min  X&Y  value 

Subroutines  Used:  MOVE 

Link:  DERIVE 


7.42  Subroutine  LSTOR 


Purpose: 

Entry  Points: 
Elements: 

Common  Areas: 
Subroutines  Used; 
Link: 


To  store  the  last  selected  points  to  the  random  file  and 
write  a flag  to  indicate  that  this  is  the  last  profile. 

None 

None 

/modelp/,Ajnits/,/idxbl/,/addres/,/data/ 

RWRITE 

SORTPR 


7.43  Subroutine  MAIN 


Purpose: 

Entry  Points: 
Elements; 
Common  Areas: 
Subroutines  Used: 
Link: 


Main  program  - controls  program  flow 

None 

None 

/GE0DET/,/0PTI0N/,AJNITS/,/PH0T0/,ADP1ND/, 

INPHDT,CORUSR,CONTR,PROFIL,SCATER,MRGE 

MAIN 
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7.44  Subroutine  MATMPY  (IAA,NRA,NCA,LAB,NRB,NCB,LAC,ICT,ICC) 


Purpose: 

Entry  Points: 
Elements: 


Common  Areas: 
Subroutines  Used: 
Link: 


To  compute  the  matrix  product  of  arrays  lAA  and  LAB. 
The  result  of  this  multiplication  is  stored  in  LAC. 

None 


I/O  Type  Vorioble  Dimension 

I R LAA  (1) 

I I NRA 

I I NCA 

I R LAB  (1) 

I I NRB 

I I NCB 

I/O  R LAC  (1) 

I I ICT 


Description 

Array  LAA 

Number  rows  in  LAA 

Number  columns  in  LAA 

Array  LAB 

Number  rows  in  LAB 

Number  columns  in  LAB 

Array  LAC  (product) 

Flag  to  designate  which  matrix 
is  to  be  transposed 


0- »R  = AB 

1- *R  = A"B 

2- *R  = AB" 

3- *R=A"  8" 


I I ICC 


Flog  to  designate  summation 
or  subtraction 


0- *R  is  stored  in  LAC 

1- *LAC+  R is  stored  in  LAC 

2- * LAC  - R is  stored  in  LAC 

3"*  - R is  stored  in  LAC 


None 

None 

MAIN 
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7.45  Subroutine  MAXM(X , Y1 , Y2) 


Purpose: 


To  test  for  slope  of  line  and  return  maximum  and  minimum  Y 
values. 


Entry  Points: 
E lements: 


Common  Areas: 
Subroutines  Used: 
Link: 


None 

I/O  Type  Vorioble 

I R X 

0 R Y1 

0 R Y2 

/CONSTA/TPDLPR/ 

XMN 

WORK 


Dimension  Description 

- X value 

- Maximum  Y value 

- Minimum  Y value 


7.46  Subroutine  MDTODD  (X,Y,X1,Y1) 


Purpose: 
Entry  Points: 
Elements: 


Common  Areas: 
Subroutines  Used: 
Link: 


To  convert  model  coordinates  to  input  digitol  data  coordinates. 
None 


I/O  Type  Variable  D imension 


/ r 

1 R 

i R 

X 

Y 

; 

Xl 

YJ 

^ Model  coordinates 

0 R 

XI 

- 

X'l 

1 Input  digital  data 

0 R 

Y1 

- 

Yj 

f system  coordinates 

Avtran/,aaram/,/Photo/ 

None 


CREATE 


Mathematics: 


Given  the  model  coordinates  of  o point,  the  computation  of  the 
digital  data  comers  involve  transformation  from  model  to  ground  and  then  ground  to 
model. 


-7  b'- 


The  transformation  from  model  to  ground  is  accomplished  by: 


X 

Y 


cosy 

-siny 


sin  y 
cosy 


Nr 

Xm 


(7.46.1) 


where 


= azimuth  relating  model  to  ground 
= model  coordinates 
= scale  relating  ground  to  model 
= origirr  of  digital  data  in  USR  coordinates 
= USR  coordinates  of  point  (>^  ,y„). 

The  digital  data  coordinates  can  now  be  computed  from  these 
computed  ground  coordinates  where 


y 

X,Y 


Xqo 

/do 


cos  O' 
sina 


-sin  a 
cos  a 


X 

Y 


1 .0/Xe^p 


(7.46.2) 


where 


a 

\/oo 
*00  /Xoo 


azimuth  relating  digital  data  to  ground 
scale  relating  ground  to  digital  data 
digital  data  coordinates  of  given  model  point. 


7,47  Subroutine  MINV  (B,N) 


Purpose: 
Entry  Points: 
Elements: 


To  compute  matrix  inversion. 

None 

I/O  Type  Vorioble  Dimension  Description 

I/O  R B (N,N)  Matrix  to  be  inverted 

I I N - Number  of  rows  and  columns  in  B 


Common  Areas: 
Subroutines  Used: 
Link: 


None 

None 

DERIVE, INPUT 
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7,46  Subroutine  MOVE  (A,B,N) 

Purpose:  Move  data  from  orroy  A to  orroy  B (does  not  destroy  data 

in  array  A) 

Entry  Points:  None 


Elements: 

I/O 

Type 

Variable 

Dimension 

Description 

I/O 

R 

A 

(1) 

Array  data  to  be  moved  from 

0 

R 

B 

(1) 

Array  data  to  be  moved  to 

1 

1 

N 

- 

Number  of  words  moved 

Common  Areas: 

None 

Subroutines  Used: 

None 

Link; 

MAIN 

7,49  Subroutine  MQFIT  (KFLAG) 

Pupose:  To  compute  multi-quad. 

Entry  Points:  None 

Elements:  I/O  Type  Voriable  Dimension  Description 

I/O  I KFLAG  - Flog  for  forming  normals 

Common  Areas;  /bATA/,/OPTION/,/STORE/,/CONST/,/f>ARMS/,/lFPDLPR/, 
Apint/,/jnits/ 

Subroutines  Used:  ADD, CLEAR, DISC, EVAL,MATMPY,MRTPLY,RREAD,RWRITE, 
SURV,TRIPLY,TRIVRT,MAXM,UNPRPK,BMAT,PROPK 
Link:  WORK 

Mathematics: 

The  solution  for  the  coefficients  representing  the  terrain  surfaces 
involves  the  solution  of  the  matrix  of  partial  derivatives  where: 

[B]  = f,^ji  (x,y).  (7.49.1) 
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7.50 


Computing  this  matrix  for  all  (x,y)  in  a block,  the  normal 
Solution  and  constant  vector  can  be  computed  from 


N = E 

and 

C = E CB]  Zj  , 

1=1 

where 

I = 1 -*  number  of  points  in  block 
Zj  = height  of  the  j th  point. 


(7.49.2) 


(7.49.3) 


The  matrix  [N]  is  formed  and  inverted  for  one  block  of  data 
only.  The  solution  of  the  coefficients  involves  the  equation 

[P]  = N'^C.  (7.49.4) 

This  [P]  is  computed  for  all  data  blocks  and  stored  to  the 

random  file. 


Subroutine  MRGE 


Purpose; 

Entry  Points: 
Elements; 
Common  Areas: 
Subroutines  Used: 
Link; 


To  combine  multiple  data  sets  for  the  final  profile  tape. 

None 

None 

AJNITS/,/t»ATA/,/PHOTO/,/MINMAX/,AANDBiy,/BUFFER/ 

MOVE,BUFFOT,RWRlTE,IMOVE,RREAD 

MERGE 


I 
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7,51  SubfooHne  MTREAD  (IWD,  KBUF,ISTAT) 


7.52 


Purpose; 


Assembly  language  subroutine  for  buffering  in  or  out  o physical 
record  from  magnetic  tope. 


Entry  Points: 

None 

Elements; 

I/O  Type 

Vorioble 

1 1 

IWD 

1 1 
0 1 

KBUF 

ISTAT 

Common  Areas: 

None 

Subroutines  Used: 

None 

Link: 

STORE,CONTR,MERC 

Dimension  Description 

- Number  of  words  to  buffer 

in  or  out 

(768)  Array  containing  data 

- Status  of  operation 


Subroutine  MRTPLY  (R,NRR,NCR,T,NRT,NCT,F,IRS) 


Purpose; 

To  mu 

Itiply 

a rectangular  matrix  and  a triangular  matrix 

Entry  Points: 

orKl  output  rectangular  matrix. 
None 

Elements; 

Type 

Variable 

Dimension 

Description 

1 

R 

R 

0) 

Rectangular  matrix 

1 

1 

NRR 

- 

Number  rows  in  R 

1 

1 

NCR 

- 

Number  columns  in  R 

1 

R 

T 

0) 

Triangular  matrix 

1 

1 

NRT 

- 

Number  rows  in  T 

1 

1 

NCT 

- 

Number  columns  in  T 

0 

R 

F 

(1) 

Resulting  matrix 

Common  Areas: 
Subroutines  Used: 
Link: 

1 1 

None 

None 

WORK 

IRS 

Option; 

-1  F = - R*T 
0F  = F+R*T 
1 F = R*T 
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7.53  SubrouHne  MULTIQ  (KNOD) 

Purpose:  To  compute  the  multiquadratic  fit  and  derive  final  profile  data. 

Entry  Points:  None 

Elements:  I/O  Type  Voriable  Dimension  Descript'cn 

I I KNOD 

Common  Areas:  /N0FIT/,/UN1TS/,/M0DELP/,/ADDRES/,/PH0T0/, 

/lASTPT/,/DATA/,/RANDBL/ 

Subroutines  Used:  CLEAR, MATMPY,MINV, LIMIT, KLEAR,MOVE,RWRlTE,RREAD 

Link:  DERIVE 

Mathematics: 

This  subroutine  will  compute  the  output  profiles  representing  the 
terrain  using  data  from  the  MMS-32  cartographic  data  tape.  On  option,  this  terrain 
will  be  represented  by  a three  to  seven  term  polynomial. 

The  coefficients  for  the  polynomial  solution  involves  the  equations: 


N = S B"  B 

i=x 

C = f B'  Z 

1=1 

where 

D = partial  derivative  of  f, 

U ~ fb  + °ay  + +04/® 

X = X,  - Xj 

y = Xi  - Xo 

Xj  ,y,  = coordinates  of  data  point  i 
Xs  lYo  - coordinates  of  local  origin 


L = N C. 


(7.53.1) 


(7.53.2) 
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The  number  of  coefficienh  usej  to  represent  o specific 
tester  area  Is  dependent  upon  the  number  of  data  points  available.  Thus,  the 
number  of  polynomial  terms  varies  from  a minimum  of  three  (3)  to  o maximum  of 
five  (5).  See  Section  3.2  for  additional  detail. 

The  fetroin  representation  is  computed  as: 

z = B'A  (7.53.3) 

where 

B = partial  derivatives 

A = computed  coefficients. 

Equation  7.53.3  is  evaluated  for  all  (x,y)  values  within  the 
area  the  coefficients  were  computed  over. 


7.54  Subroutine  MVSTR  (IL,N'BLTP,FELEV) 


Purpose: 
Entry  Points: 
Elements: 


Comnxsn  Areas: 
Siibroutines  Used: 
Link: 


To  locate  and  store  data  with  respect  to  specific  profiles. 

None 

l/Q  Type  Variable  Dimension  Description 

^ ^ IL  - Data  index  within  physical  block 

I I NBLTP  - Physical  data  block  count 

^ ^ FELEV  - Elevation  associated  with  current 

data 

/modelp/,/index/,Ajnits/,/Photo/,/data/,/lvtran/ 

STOR 

SORTPR 
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Mathematics: 

Having  read  a block  of  data  from  the  MMS-32  contour  tape,  oil 
coordinates  must  be  converted  to  the  model  coordinate  system  for  comparison  to 
desired  profile  coordinates.  If  these  coordinates  fall  within  a boundary  surrounding 
a profile,  the  values  ore  stored. 

The  transformotion  from  MMS-32  coordinotes  to  model  coordinates 


Involve: 


— 

- — 

- — 

X 

- 

cos  (^Z) 

sin  (AZ) 

Xoo 

1 

0 

0 

Xq  «a 

Y 

— _ 

-sin  (AZ) 

cos  (AZ) 

•< 
0 
^ 0 

Yq  BO 

(7.54.1) 


where 


AZ 

^/oo 
^«8  /Xj  R 0 


azimuth  relating  digital  data  to  ground 
digital  data  coordinates  of  point 
scale  relating  ground  t digital  data 
USR  coordinates  of  digital  data  origin. 


where 


cos(4.ZlM)  -sin  (^ZIM) 
sin(^ZIM)  cos(^ZIM) 


• 1000/x«^^  ! 


AZIM 

X,Y 


azimuth  relating  ground  to  model 
USR  coordinates  of  point 
scale  relating  ground  to  model. 


7.55  Subroutine  MXMN  (X,YMX,YMN) 


I 

7.56 


Puqiose: 
Entry  Points: 
Elements: 


Common  Areas: 
Subroutines  Used: 
Link: 


To  compute  minimum  and  maximum  Y for  each  block. 
None 

I/O  Type  Vorioble  Dimension  Description 

I R X - X input 

0 R YMX  - Maximum  Y 

0 R YMN  - Minimum  Y 

/Tpdlpr/,Ajnits/ 

XMN 

PROFIL 


Subroutine  NEW 


Purpose:  To  compote  the  new  profiles. 

Entry  Points:  None 

Elements:  None 


Common  Areas: 


Subroutines  Used: 
Link: 


/PARMS/,AMNMAX/,/FPDLPR/,AHOTO/,/t)ATA/, 

Ainits/,/const 

DMATX,MATMPY,PPACK, CLEAR, COFF 
EXTRA 


AAathematics: 

The  computed  height  of  the  terrain  at  point  (x,y)  is  computed  by 

= f7.ia(X/y)  (7.V 

Z = [B]  • [P] 

where 

CB]  = matrix  of  partial  derivotives  evaluated  fa 
[P]  *=  computed  coefficients  representing  t<- 
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/'AD-A043  456 


OFF-LINE  ORTHOPRINTER  SOFTWARE(U)  DBA  SYSTEMS  INC 
MELBOURNE  FLA  d MAUPIN  ET  AL.  AUG  77  RADC-TR-77-256 
F30602-74-C-0297 


UNCLASSIFIED 


F/G  9/2 


NL 


1 


7.57  Subroutine  OAMTIV  (ANG,OMATRX) 


Purpose:  To  compute  L.V.  orientation  matrix  given  L.V.  angles. 

Entry  Points;  None 

Elements:  I/O  Type  Vorioble  Dimension  Description 

I R ANG  (3)  L.V.  angles  a 
0 R OMATRX  (9)  Arroy  containing  computed 

orientation  matrix 

G>mmon  Areas;  None 

Subroutines  Used:  None 

Link:  SDT 

Mathematics: 

The  ground  to  photograph  local  vertical  orientation  matrix  is 

defined  as: 

Ttr 

where 

T|.v(lrl)  “ *<n  Of  cosfp 
= cosasinfp 
Tu,0,3)  = sin<ip 

T|,y(2/1)  = cos  a cos  u- sin  a sin  <p  sin  0) 

T^y(2,2)  =*ein  acosbj- cosofsin  <psin  0)  (7.57.1) 

T^y(2,3)  = cosfpsinw 

T|.y(3,l)  = -cos a sin  0)  - sin  a sin  (p  cos  w 
T^y(3,2)  s sin  a sin  CO- cos  a sin  ^ cos  10 
T^y(3,3)  = cosipeosu) 
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7.58  Subroutine  OUTPRO  (K,NOY) 


Purpose: 

Outputs  profile  data  to  file  colled  tOUT . 

Entry  Points: 

None 

Elements: 

I/O  Type  Variable  Dimension  Description 

II  K - Profile  number 

1 1 NOY  ” Total  number  of  Y coordinate 

values  in  profile 

Common  Areas: 

/lDXBl/,/MODELP/,/UNITS/,/DATA/ 

Subroutines  Used: 

IMOVE 

Link: 

SORTER 

7,59  Subroutine  PACK  (K,FX,FY,FZ,IDX) 


Purpose: 
Entry  Points: 
Elements: 


Common  Areas: 
Subroutines  Used: 
Link: 


Fo  move  o 32-bit  real  word  into  two  16-bit  integer  words. 
None 


'bATA/ 

None 

:ONTR 


[ffi® 

Voriable 

Dimension 

Description 

1 

K 

Profile  ID 

R 

FX 

- 

X value 

R 

FY 

Y value 

R 

FZ 

• 

Z value 

R 

IDX 

- 

First  address  of  point  in 

common  block 
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7.60  Subroutine  PK  (VALE,IFX) 


Purpose:  To  pack  one  real  word  into  two  integer  words. 

Entry  Points:  None 

Elements:  I/O  Type  Vorioble  Dimension  Description 

I R VALE  - Real  word  input 

0 I IPX  (2)  Integer  word  array 


Common  Areas:  None 

Subroutines  Used:  None 

Link:  STORE 


7.61  Subroutine  PLHXYZ  (PLH,XYZ) 

Purpose:  To  convert  geographic  coordinates  of  geocentric  coordinates 


Entry  Points: 

None 

Elements: 

I/O  Type 

Variable 

Dimension 

Description 

1 R 

PW 

(3) 

Geographic  coordinates 

Common  Areas; 

0 R 

/GEODET/ 

XYZ 

(3) 

Geocentric  coordinates 

Subroutines  Used:  None 

Link:  INPUT 
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7.62  Subroutine  PPACK  (1,X,Y,Z,K) 


7.63 


7.64 


i 


Purpose: 


To  pock  real  words  in  two  interger  words. 


Entry  Points: 

None 

Elements: 

Type 

Variable 

1 

1 

1 

1 

R 

X 

1 

R 

Y 

1 

R 

Z 

1 

1 

K 

Dimension  Description 

- First  location  of  real  words 

In  integer  array 

- X value  to  be  stored 

- Y value  to  be  stored 

- Z value  to  be  stored 

- K offset  from  first  location 

in  integer  array 


Common  Areas:  /DATA/ 

Subroutines  Used:  None 


Link:  EXTRA 


Subroutine  PROFIL 


Purpose: 
Entry  Points: 
Elements: 


Main  program  for  processing  profile  data. 

None 

None 


Common  Areas:  /0PTI0N/,/C0NST/,/PARMS/,/UNITS/,/FPDLPR/, 

/Store/,Apint/,/minmax/,/nodes/,/photoa 

/PARAM/,/LVTRAN/ 


Subroutines  Used:  BOUND, FPB0ND,GEN0DE,IRANI0,MQFIT,PR0FLE,ASTRP 
Link:  PROFIL 


Subroutine  PROFLE 


Purpose; 

Entry  Points: 
Elements: 

Common  Areas: 

Subsfeutines  Used: 
Link: 


To  compute  final  output  profiles. 

None 

None 

/UNITS/,/t>ATA/,/MlNMAX/,/FPDLPR/,/1>HCTO/, 

/const/,/param/,Avtran/ 

CONTF,PROPAK,XYINT 

CREATE 
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I 


7.65  SubfooMne  PROJ  (XY,PHLAHT) 


Purpose:  To  colculote  universal  space  rectangular  coordinates 

of  the  point  of  intersection  of  a single  ray  from  a 
spatial  photograph  with  a reference  ellipsoid. 

Entry  Points:  None 

Elements:  I/O  Type  Vorioble  Dimension  Des' 

I R XY  (2)  Image  coor 

0 R PHIAHT  (3)  of 

Common  Areas:  /GE0DET/,/0PTI0N/,/PARAM/,/UNPAR/ 

Subroutines  Used:  MATMPY,MOVE,OMATLV,XYZPLH,USRLVR 


I/O  Type  Vorioble  Dimension 

1 R XY  (2) 

0 R PHIAHT  (3) 


Description 

Image  coordinate  of  point 
<p,X,H  of  imoge  input 


Mathematics: 

Given  the  image  coordinates  and  focal  length  of  the  camera, 
the  direction  cosines  can  be  computed  by; 


V = Tu 


where 

x,-y  *=  image  coordinates  of  point 
-f  = focal  length  of  camera 
Ty,j  = USR  orientation  matrix. 

Knowing  the  direction  cosines,  the  USR  coordinates  can  be  computed  by; 

Z(l)  = (-B+D)/2A 
Z(2)  = (-B-D)/2A 

X(l)  = UAV(Z1)+ (X  + X(t))-U/W(t) 

X(2)  = U/W(Z2)  + (X+X(t))-U/W(t) 

Y(l)  = V/W(Z1)+  (Y+Y(t))-V/W(t) 

YC2)  = VA/  (Z2)  4 (Y4  Y(t))  - V/V/(t) 
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wher* 


A » (U/W)P  + + 1/(1 " 

B - 2lM^ ((X+  X(t))-  U/W(Z  + Z(t)))  + (V/W)((Y+  Y(t))  - V/Y(Z - Z(t)))] 

C = [X  + X(t)  - U/W(Z  - Z(t))?  + [Y+  Y(t)  - VA»/(Z  - Z(t))]=  - a» 

i 

D = [B»  -4AC] 

e*  * eccentricity 

t •=  time  from  center  of  scon  (For  pan) 

X,Y,Z  s exposure  station  coordinotes 

• e e 

X,Y,Z  = linear  rates 
a “ semi-major  axis. 


From  the  above  computed  X,Y,Z  values,  the  distance  from  the  origin  can  be 


computed  by: 

R(l)  *=  C(X(l)-X:)8+(Y(l)  -Yo)a4(Z(l)-Zo)*]^ 
R(2)  - [(X(2)  - Xo  )*  + (Y(2)  - Yo)»  + (Z(2)  - Zo)= 


(7.65.3) 


where 

/Yg  ,Zg  = coordinates  of  the  exposure  station . 


The  comparison  of  the  R values  determine  desired  point  of  intersection.  The  USR 
coordinates  will  be  X(1),Y(1),Z(1)  unless  R(1)  ^R(2),  then  the  USR  coordinates  X(2),Y(2),Z(2) 

ore  used. 


-89- 


7.66  Subroutine  PROPAK  (K,X,Y,Z) 


Purpose: 

To  pock  o real 

word  into  common  area  occupied  by  two  integer 

words. 

Entry  Points; 

None 

Elements: 

I/O  Type  Vorioble  Dimension 

Description 

1 1 

K 

Point  location  in  array 

1 R 

X 

X value 

1 R 

Y 

Y value 

1 R 

Z 

Z value 

Common  Areas: 

/Data/ 

Subroutines  Used: 

None 

Link: 

CREATE 

Subroutine  PROPK 

(|,X,Y,Z,K) 

Purpose:  To  pock  o real  word  into  oreo  occupied  by  two  integer  words. 


Entry  Points:  None 


Elements: 

I/O  Type 

1 1 

1 R 

1 R 

1 R 

1 1 

Variable  Dimension 

1 

X 

Y 

Z 

K 

Description 

Indicates  point  location  in  orray 
X value  to  store 
Y value  to  store 
Z value  to  store 
Offset  of  storage  location  from 
first  location  in  array 

Common  Areas: 

A>ATA/ 

Subroutines  Used: 

None 

Link: 

PROFIL 
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7.68  Subfootine  PTAPE 
Purpose: 


Prepare  data  for  punching  AS  1 1 shutdown  tope . 


Entry  Points:  None 

Elements:  None 


Common  Areas:  AJN1TS/,/0PTI0N/,AJNPAR/,/GE0DET/,/GR0UND/, 

/PARAM/,/PHOTO/,/OFF/ 

Subroutines  Used:  ARATES,ASllAN,CLEAR,MATMPY,MOVE,OMATLV, 

PLHXYZ,PUNDEC,USRLVR,PUNTAP 

Unk:  SDT 


Mothemotics: 


where 


Compute  Y axis  in  LV  system 


r..r  n 

Xg 

Xj  - 

Y» 

Yx-Ye 

3. 

Xo  ,Yo  ,Zo 
X,,Y,,Z, 


ore  USR  values  of  projected  Y axis 
USR  values  of  origin 
LV  values  of  projected  Y axis 

**  ^7.94  • 


(7.68.1) 


Compute  azimuth  of  model  system 
AZIM  = ore  ton  [Xg/V,  ] . 

Compute  By  ,6^  ,B, 


“b." 

"x. 

Br 

K P 

’'uss  LSS 

Y.  -Y, 

Ll* 

Z«-Zo 

(7.68.2) 


(7.68.3) 
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where 


/Y|  ,Z(  - coordinates  of  exposure  station 
r®r  /®i  ~ components. 


The  otmospheric  refraction  coefficients  are  computed  using 
a polynomiol  derived  to  fit  the  graphic  table  shown  in  Figure  6-1  and  Figure  6-2 
of  the  volume  I of  the  AS-1  lB-1  manuals.  If  the  flying  height  is  above  50  kilometers 
a maximum  set  of  values  are  used.  For  under  50  kilometers  the  following  equation  is 
used  to  compute  the  AS-1 1 atmospheric  coefficients  (CE  and  CH). 


CH  = (q,  +aiH+03H»  +a3H3+o^H«  + agHS)(l/MS)(C) 
CE  = Ofe  +bjH +bjH=  + bgH® +b^H*  + bBH®)(MS/1000) 

where 


'MS  = Model  Scale 

H = flying  Height  of  vehicle  (in  kilometers) 


c 

S= 

l.OE  + 06 

s 

-243.001 

lb 

= 1.698E-03 

s: 

115.042 

bs 

= -3.943E-04 

s 

- 7.168 

bs 

= 4.030E-05 

Oj 

=: 

18.5999 

ba 

= -1.911E-06 

ss 

- .8107 

b* 

= 4.067 

<% 

cs 

.00928 

bs 

= -3.133 

(7.68.4) 


7.69 


Subroutine  PUNCH  (ICH) 

Purpose: 

Assembly  language  routine  for  punching  paper  tapes. 

Entry  Points: 

None 

Elements: 

I/O  Type  Variable  Dimension  Description 

1 1 ICH  - Character  to  be  punched 

Common  Areas; 

None 

Subroutines  Used: 

None 

link: 

SDT 

-92- 


7.70  Subrootine  PUNDEC  (DNUM,NODEC) 


Purpose: 

To  punch  dota  on  paper  tape. 

Entry  Points: 

None 

Elements: 

I/O  Type  Varioble  Dimension 

Description 

I R DNUM 

Decimal  number  to  be  punched 

1 1 NODEC 

Number  of  decimal  places 

Common  Areas: 

/UNITS/ 

Subroutines  Used: 

CLEAR 

Link: 

SDT 

Subroutine  PUNTAP 

Purpose: 

Paper  tape  punch  routine. 

Entry  Points; 

None 

Elements: 

None 

Common  Areas: 

Ajnits/ 

Subroutines  Used: 

PUNCH 

Link: 

SDT 

Subroutines  PXYZ  (PLH,XYZ) 

Purpose: 

Transform  geographic  coordinates  to  geocentric  coordinates. 

Entry  Points: 

None 

Elements: 

I/O  Type  Variable  Dimension 

Description 

1 R Pm  (3) 

Geographic  latitude, longitude 
and  height 

0 R XYZ  (3) 

USR  coordinates 

Common  Areas: 

/GEODET/ 

Subroutines  Used: 

None 

Link: 

EXTRA 
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7.73  Subroutine  RASH 


Purpose: 

Entry  Points: 
Elements: 

Common  Areas: 
Subroutines  Used: 

Link: 


Input  AS  11  shutdown  tope  data  and  compute  rotation  matrix. 

None 

None 

/param/,AJNits/,/geodet/ 

AS  1 1M,UTMPRP,UTMGE,USRLVR,  PLHXYZ,MATMPY,XYZPLH , 
MOVE 

CORD 


7. 74  Subroutine  RGEO 


Purpose: 

Entry  Points: 
Elements: 
Common  Areas: 
Subroutine  Used: 
Link: 


To  create  additional  point  file.  Data  read  in  geographies. 

None 

None 

A)  NITS/ 

DMSRAD,PLHXYZ 

SDT 


Subroutine  RREAD  (IUN,IFWA, 

NWDS,IDIS) 

Purpose: 

To  read  from  the  random  file. 

Entry  Points: 

None 

Elements: 

I/O  Type 

Variable  Dimension 

Description 

I/O  1 

lUN 

Unit  reading  from 

1 1 

IFWA 

First  location  to  read 

1 1 

NWDS 

Number  of  words 

1 1 

IDIS 

Storage  location  on  random  file 

Common  Areas: 

AlAND  Bl/,/UN  ITS/,/DATA/ 

Subroutines  Used: 

None 

Unk: 

PROFIL,EXTRA, MERGE  ,CONTR 

( 
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7.76  Subroutine  RTPLY  (R,NRR,NCR,T,NRT,NCT,F,IRS) 


7.77 


Purpose: 

To  multiply  a rectangular  matrix  and  a triangular  matrix. 
Output  is  a rectangular  matrix. 

Entry  Points: 

None 

Elements: 

I/O  Type  Variable  Dimension  Description 

1 R R (1)  Rectangular  matrix 

1 1 NRR  - Number  rows  in  R 

1 1 NCR  - Number  columns  in  R 

1 R T (1)  Triangular  motrix 

1 1 NRT  - Number  rows  in  T 

1 1 NCT  - Number  columns  in  T 

0 R F (1)  Resulting  matrix 

1 1 IRS  - Optiort.: 

-I  -R*T 

0 C+R*T 

1 R*T 

Common  Areas: 

None 

Subroutines  Used: 

None 

Link: 

EXTRA 

Subroutine  RUSR 

Purpose: 

To  create  additional  point  file.  Data  read  in  USR 
coordinates. 

Entry  Points: 

None 

Elements: 

None 

Common  Areas: 

AJNITS/ 

Subroutines  Used: 

None 

Link: 

SDT 
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7.78 


Subroutine  RUTM  (OPHI,OLAM) 


Purpose: 

Entry  Points: 
Elements: 


Common  Areas: 
Subroutines  Used: 
Link: 


To  read  in  additional  points  in  UTM  coordinates  and  output 
on  a nie  their  USR  coordinates. 

None 

Dimension  Description 

Origin  of  the  UTM  values 

Ajnits/,/T>hoto/,/geodet/,/option/ 

UTMPRP,UTMGE,PLHXYZ 

CORD 


I/O  Type  Vorioble 

I R OPHI 
I R OlAM 


7.79  Subroutine  RWRITE  (IUN,IFWA,NWDS,IDIS) 


Purpose: 

To  write  to  the  random  file 

Entry  Points: 

None 

Elements: 

I/O  Type.  Variable  Dimension  Description 

I/O  1 lUN  - Unit  to  write  on 

1 1 IPA/A  - First  location  in  common 

1 1 NWDS  - Number  of  words 

0 1 IDIS  - Storage  location  on  random  file 

Common  Areas: 

AJN  ITS/,/RA  ND  Biy,/DATA/ 

Subroutines  Used: 

None 

Link: 

PROFIL,EXTRA,MERGE,CORTR 

Subroutine  SCATER 

Purpose: 

Control  program  flow  when  using  scattered  points  for  making 
an  orthophoto. 

Entry  Points: 

None 

Elements: 

None 

Common  Areas: 

AJN1TS/,/GE0DET/,/DATA/,APDLPR/,/N0DES/,/C0NST/, 

/Photo/,/Param/,/minaaax/,/parms/,/Tpint/ 

Link: 

EXTRA 
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7.81 


Subroutine  SOLVE 


Purpose: 

Entry  Points: 
Elements: 

Common  Areas: 

Subroutines  Used: 
Link: 

MothemoticS: 


Form  and  solve  the  normal  equations  and  compute  the 
coefficients. 

None 

None 

/NODES/,/t)ATA/,/CONST/,AARMS/,/FPDLPR/,A<NITS/, 
API  NT/ 

BMATX,TRIPLY,MATMPY,TIVRT,RTPLY,CLEAR 

EXTRA 


The  solution  for  the  coefficients  representing  the  terrain  surface 
Involves  the  solution  of  the  matrix  of  portiol  derivatives  were 


CB]  = f7.i3  (x,y). 


(7.81.1) 


Computing  this  matrix  for  all  (x,y),  the  normal  solution 
and  constant  vector  can  be  computed  from 


N = E [B][B^]  (7.81.2) 

ond 

C * E [B]  Z (7.81.3) 

where 

J *=  1 -*  number  of  points 
Zj  *=  height  of  the  ) th  point. 

The  matrix  [N]  is  found  and  inverted  for  one  block  of  data 
only.  The  solution  for  the  coefficients  involves  the  equation 

[P]  = N“'c.  (7.81.4) 
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7.82  Subroutine  SORTBL 


Purpose; 

Entry  Points: 
Elements; 

Common  Areas: 
Subroutin  fs  Used: 
Link: 


To  read  data  from  random  file  and  sort  as  to  profiles. 

None 

None 

/lDXBiy,AJNITS/,/ADDRES/,/MODELP/,/ADIDBl/,/DATA/ 

KLEAR,OUTPRO, PACK, RREAD,RWRITE, UNPACK 

SORTPR 


7.83  Subroutine  SORTPR 


Purpose: 

Entry  Points: 
Elements: 

ComiTK>n  Areas: 

Subroutines  Used: 
Link: 


To  read  in  sets  of  four  profiles. 

None 

None 

AiNITS/,/LASTPT/,/NOFIT/,/MODELP/,/DATA/,/OPTION/, 
/^NDBL/  ' 

KLEAR,INCLUD,SRFILL 

DERIVE 


7,84  Subroutine  SRFIL 


Purpose: 

Entry  Points: 
Elements: 

Common  Areas: 
Subroutines  Used: 
Link: 


To  sort  profile  sets  into  raster  data. 

None 

None 

Aastpt/,/Nofit/,AJNits/,/modelp/,/data/ 

MULTIQ 

DERIVE 
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7.85  Subroutine  STOR  (K,FX,FY,FZ) 


Purpose: 
Entry  Poinis: 
Elements: 


Common  Areas: 
Subroutines  Used: 
Link: 


To  store  selected  points  of  contour  data  to  a random  file  . 
None 


Im 

Variable 

Dimension 

Description 

1 

1 

K 

Profile  ID 

1 

R 

FX 

- 

Profile  X value 

1 

R 

FY 

- 

Profile  Y value 

I 

R 

FZ 

- 

Profile  Z value 

/mode  LP/,  AJ  N ITS/,/IDX  Bl/,/ADDRES/,/bATA/ 

PACK,RWRITE 

SORTPR 


7.86  Subroutine  SURV  (C) 

Purpose:  To  include  random  point  contribution  in  normal  equation. 


Entry  Points: 
Elements: 


None 


I/O  Type  Variable  Dimension  Description 
I/O  R C (70)  Constant  vector 

Common  Areas:  /C0NST/,APDLPR/,/DATA/,/UNITS/ 

Subroutines  Used:  BMAT,MATMPY,TR1PLY 

Link:  WORK 

Mathematics: 

Given  image  coordinates  of  the  randomly  distributed  survey  or 
terrain  points,  the  partial  derivative  matrix  is  computed  os: 


(7.86.1) 


This  [B]  matrix  contribution  is  summed  into  the  normal  equation 
and  constant  vector  formed  in  the  profile  block  (Section  7.50). 


I 


I 
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N 


(7.86.2) 


•=  E [B3[Bn 

C »=  L [B]Zj  (7.86.3) 

where 

J *s  1 -»  total  number  of  points. 

These  matrices  ore  used  in  the  subroutine  described  in 
Section  7.49.  They  ore  entered  os  the  solution  to  equations  7.49.2  and  7.49.3. 


7.87  Subroutine  TIVRT  (TMT,SVC,MAR,MFI) 


Purpose: 


Inversion  and/or  solution  for  a triangular  matrix. 


Entry  Points; 
Elements; 


Common  Areas: 
Subroutines  Used: 
Link: 


None 


1^ 

IZEi 

Variable 

Dimension 

Description 

I/O 

R 

TMT 

(1) 

Matrix  stored  upper  triangular 

I/O 

R 

SVC 

(1) 

Solution  vector 

1 

1 

MAR 

- 

Matrix  rank 

I 

I 

MFi 

Options: 

-1  inverse  only,  no  solution 

0 solution  only,  no  inverse 

1 both  inverse  and  solution 

None 

None 

EXTRA 
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7.88  Subroutine  TR  (X,Y,Z,Z1,XX,YY,ZZ) 


Purpose:  To  transform  points  from  USR  coordinates  to  model  scale. 


Entry  Points:  None 


Elements: 

hEl 

Variable 

Dimension 

Description 

1 

R 

X 

- 

X 

USR  coordinate 

1 

R 

Y 

- 

Y 

USR  coordinate 

1 

R 

Z 

- 

Z 

USR  coordinate 

1 

R 

Z1 

- 

z 

UTM  value 

0 

R 

XX 

- 

X 

model  coordinate 

0 

R 

YY 

- 

y 

model  coordinate 

0 

R 

ZZ 

- 

z 

model  coordinate 

Common  Areas: 

/Photo/,/param/ 

Subroutines  Used: 

PXYZ,MATMPY 

Link:  EXTRA 

Mothematics: 

Given  the  geographic  coordinates  of  the  origin,  the  USR 
coordinates  can  be  computed  from: 


(7.88.1) 


Knowing  the  USR  coordinates  of  the  origin  and  the  USR 
coordinates  of  a point,  the  model  coordinates  can  be  computed  by: 


"Xg  - X," 

Yu 

Yu  -Yo 

_Zt 

_Zg  -Zo_ 

where 

Xu,Yg,Zy  = USR ':oordinates  of  point 
^s*  LIS  ~ rotation  from  USR  to  LSR  system. 


(7.88.2) 
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(7.88.3) 


where 


cos  (AZ)  sin  (AZ)  0 

-sin  (AZ)  cos(^Z)  0 

Yl 

0 0 1 

_Zt 

AZ 


azimuth  relating  ground  to  model 
scale  relating  ground  to  model. 


7,69  Subroutine  TRANFR 


Purpose: 

Entry  Points: 
Elements: 

Common  Areas: 
Subroutines  Used: 
Link: 


To  transform  data  from  USR  coordinates  to  model  coordinates. 

None 

None 

Ahoto/,Ainits/,/const/,/param/ 

MOVE,PXYZ,MATMPY,TR 

EXTRA 


7.90  Subroutine  TRANS  (XY,XYZUSR) 


Purpose: 
Entry  Points: 
Elements: 


Common  Areas: 
Subroutines  Used: 
Link: 


Compute  a two-dimensional  transformation. 

None 

I/O  Type  Voriable  Dimension  Description 

I R XY  (2,4)  Image  coordinate  of  D.D.  corners 

I R XYZUSR  (3,4)  Object  coordinate  of  D.D.  corners 

Ainits/,/geodet/,/photo/,Autran/ 

CLEAR  ,MATMPY,MINV, MOVE,  PLHXYZ,USRLVR 
COR 
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I 

» 

J 


Mathematics: 

Given  the  equations 

f,  *=  ( xcostt+y  sin  a)S -X,  (7.90.1) 

f,  *=  (-X  sin  a+  y cosa)S  - (7.90.2) 

where 

x,y  = Image  coordinates 

a *=  azimuth 

S B scale 

Xj » Yp  = LSR  coordinates  of  comers . 


To  solve  equotions  7.90.1  and  7.90.2  on  iteration  of  the 
least  square  solution  is  performed . 


The  matrix  of  partial  derivatives  of  fj^  ond  f^  is  computed. 
The  elements  comprising  the  partial  derivative  matrix  is: 


CB] 


cc] 


6f,/ba  6f,/JXp  6f,/aYp  6f,/^S 

Bf,/6a  6f/aXp  bf/5Yp  6f/aS 

The  constant  vector  is  comprised  of 


X - ( X cosa+  y sin  a)S  +Xp 
Y - (-X  sin  a + y costt)S  + Yp 


# 


(7.90.3) 


(7.90.4) 


The  least  squares  solution  involves 

4 

N *=  S 
1=1 

C = E B^C 
1=1 

A = N“'c 


(7.90.5) 

(7.90.6) 

(7.90.7) 
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rtie  values  contained  in  the  [A]  matrix  ore  corrections  for  a,Xj  ,S  where 


a„  B Oo  +£Aof 
X„  - X,o+EAX„ 

Y„  = Y,„  + LAY, 

S,  = Sb  +EAS, 

where 

Aa,AXg  ,AYj  ,AS  = change  in  a ,Xj, , Y^  ,S 
«o  *^eo  r'*^oo  “ initial  approximations 
0C„  «^SN  '^ON  ~ updated  values. 


(7.90.8) 


7.91  Subroutine  TRIPLY  (LAA,NRA,NCA,LAB,NRB,NCB,LAC,ICT,ICC) 


Purpose: 


Entry  Points: 


To  compute  the  matrix  product  of  arrays  LAA  and  LAB 
storing  the  upper  triangular  portion  only.  The  result 
of  this  multiplication  is  stored  in  LAC. 

Norte 


Elements: 

Voriable 

Dtmertsion 

Description 

1 

R 

LAA 

0) 

Arroy  LAA 

1 

1 

NRA 

- 

Number  rows  in  LAA 

1 

1 

NCA 

- 

Number  columns  in  LAA 

1 

R 

LAB 

(1) 

Array  LAB 

1 

1 

NRB 

- 

Number  rows  in  LAB 

1 

1 

NCB 

- 

Number  columns  in  LAB 

I/O 

R 

LAC 

(1) 

Arroy  LAC  (product) 

1 

1 

ICT 

FLAG  to  designate  which 
matrix  are  to  be  transposed: 

0- *R  = A B 

1- ♦R  = A"B 

2- *R  = A B" 

3.»r  = A"B' 

I I ICC 


Common  Areas:  None 

Subroutines  Used:  None 

Link:  MAIN 


FLAG  to  designate  summation 
or  subtraction 

6 -tR  is  stored  in  LAC 

1 •♦LAC  + R is  stored  in  LAC 

2 -♦LAC  - R is  stored  in  LAC 

3 •♦  - R is  stored  in  LAC 
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7.92  Subroutine  TRIVRT  CTMT,SVC,MAR,MFI) 


Purpose: 

Inversion  solution  for  o 

triangular  matrix. 

Entry  Points: 

None 

Elements: 

I/O  Type 

Variable 

Dimension 

Description 

I/O  R 

TMT 

0) 

Matrix  stored  upper  triangular 

I/O  R 

SVC 

(1) 

Solution  vector 

1 1 

MAR 

- 

Matrix  rank 

Common  Areas: 
Subroutines  Used: 
Link: 

1 1 

None 

None 

WORK 

MFI 

FLAG  to  designate  solution 
and/or  inversion 
-1  - inverse  only,  no  solution 

0 - solution  only,  no  inverse 

1 - both  invert  + solution 

7.92  Subroutine  UNPACK  (K,FX,FY,FZ,1DX) 


Purpose: 


To  move  two  16-bit  integer  words  to  a 32-bit  reel  word* 


Entry  Points: 
Elements: 


Common  Areas: 
Subroutines  Used: 
Link: 


None 


J/O 

Type 

Variable 

Dimension 

Description 

0 

1 

K 

- 

ProHle  ID 

0 

R 

FX 

- 

X value 

0 

R 

FY 

- 

Y value 

0 

R 

FZ 

- 

Z value 

1 

1 

IDX 

- 

Address  of  point  in  common  block 

/DATA/ 

None 

CONTR 
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7.94 


7.95 


Subroutine  UNPRPK  (K,L,X,Y,Z) 

Purpose: 

To  unpock  real  words  which  have  been  stored  in  two 
integer  word  locations. 

Entry  Points: 

None 

Elements: 

I/O  Type  Variable  Dimension  Description 

I 1 K - Word  location  in  common  > 

II  L - Number  to  indicate  words 

OR  X - X value 

OR  Y - Y value 

OR  Z - Z value 

Common  Areas: 

/bATA/ 

Subroutines  Used: 

None 

Link: 

PROFIL 

Subroutine  USRLVR 

(PLG,ROTMAT) 

Purpose: 

To  compute  transformation  from  geocentric  to  local  space 
rectangular. 

Entry  Points: 

None 

Elements: 

I/O  Type  Variable  Dimension  Description 

1 R PLG  (3)  Geographic  <p,X,y of  LSR 

0 R ROTMAT  (3,3)  Transformation  matrix 

Common  Areas: 

None 

Subroutines  Used: 

None 

Link: 

INPUT 

Mothematics: 

Gtven  the  geographic  coordinates  for  the  origin  of  the  LSR 
system  and  the  azimuth  of  the  "Y"  oxis  from  north  the  rotation  matrix  from  USR  to 
local  space  rectangular  (LSR)  is  defined  as:  (7.95.  T) 


>co$ys!n  X + sin  V sin  tp  ccaX 
cosy  cosX  -t  sin  y sin  <psin  X 
- sin  y cos<p 


-sin  y sin  X - cosy  sin  cos X cos«pcosX 

sin  y cosX  - cosy  sin  ^ sin  X cos<psinX 

cosycossP  sin^ 
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7.96  Subroutine  UTMGE  (OPHI  ,OlAM ,XG ,YG ,OK,OUT,PHI  ,ALAM) 


Purpose: 

Tratuform  UTM  to  geodetic 

coordinates. 

Entry  Points; 

None 

.Elements: 

]/0  T^ 

Vorioble  Dimension 

Description 

1 R 

OPHI 

- 

UTM  origin  on  map  sheet  (c) 

1 R 

OLAM 

- 

UTM  origin  on  mop  sheet  (X) 

1 R 

XG 

- 

UTM  coordinotes  to  be  transformed 

1 R 

YG 

- 

UTM  coordinates  to  be  tronsformed 

1 R 

OK 

- 

Scale  factor 

1 R 

OUT 

(9) 

Matrix  necessory  for  UTM  conversion 

0 R 

PHI 

- 

^of  converted  coordinate 

Common  Areas: 
Subroutines  Used: 
.Link: 

Mathematics: 

0 R 

/GEODET/ 

None 

CORD 

AlAM 

X of  converted  coordinate 

Conversion  process  for  converting  Universal  Transverse 
^creator  to  Geographies 


X *=  0«, -BJ/OK  (7.96.1) 

y = Y./OK 

-^rbere 

X,  ,Yj  = UTM  coordinates  of  point  G 
OK  = scale  factor 

*=  air  base  component  in  X 

R = o(l  - ^)/(l  - c®  sin®  <Po (7.96.2) 
■where 

o « semi-major  axis 
e?  * eccentricity 

^ *5  UTM  origin  in  <p  of  Geographies 

9>‘-  «Po  + yA  (7.96.3) 
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A = OLJT{3) 

B = OUT(4) 

C = OUT(5)  (7.96.4) 

D «=  CUT(6) 

where 

OUT(9)  = nine  element  matrix  from 
UTMPRP  Section  7.96) 


S = a(l-e®)[A(<p'-<PQ)-B(sin2<p*-sin2<Po)  + (7.96.5) 

C($in4^'-sin4<(55)  - D(sin6tp‘-sin6<pQ)] 

«p'=  <p'-(s-y)A  (7.96.6) 

Return  to  equation  7.96.5  until  s-y  < .005  or  until  ten  iterations  are  made. 

N = a/(l-e®  $in=<p')^  (7.96.7) 

t?“=  cos*<p'[e®/l -e^]  (7.96.8) 

R = a(l -e®)/(l -e®  sin®  <p')®/^  (7.96.9) 

tanip' C-x?/2R(X/KI)  + x*/24RN®(5+3  tan®<p')]  (7.96.10) 

tp  B ip*  + ^ (7.96.11) 

AX  = l/cos<p'[X/M  - l/6(XAl)®(l+2tan®<p'+Tj®)  (7.96.12) 

+ 1/120  (V^J)®  (5  + 28  tan®  ^*  + 24  tan*  (p')] 

X = Aq  - ax  (X  positive  west  of  Greenwich).  (7.96.13) 


7.97  Subroutine  UTMPRP  (OPHI,OUT) 


Purpose: 

Entry  Points: 
Elements: 

Common  Areas: 
Subroutines  Used: 
Link: 

Mathematics: 
transforming  UTM 

OUTO)  = 
OUT  (2)  = 
OUT  (3)  = 
OUT  (4)  = 
OUT(5)  = 
OUT(6)  = 
0UT(7)  = 
OUT(8)  = 
OUT(9)  = 


Data  preparation  for  computing  USR  coordinates  from 
UTM  coordinates. 

None 

I/O  Type  Vorioble  Dimension  Description 

I R OPHI  - PHI  of  origin  of  UTM  coordinates 

0 R OUT  (9)  Matrix  necessory  for  UTM  conversion 

/GEODET/ 

None 

CORD 

This  routine  forms  a nine  element  array  of  constants  needed  for 
coordinates  to  geographies 

sin  <p 
cos(p 

^ + (3/4)  e®  + (45/64)  e*  + (175/256)  e® 

1/2  [(3/4)  e®  + (15/16)e“  + (525/512)^  J (7.97.1) 

1/4  [(15/64) e*  + (105/256) e®] 

1/6  [(35/512)^] 
sin  2cp 
sin  4(fi 
sin  6p. 


7.98  Subroutine  XFILL  (IPR,Z,X,J) 


Purpose:  To  fill  block  in  X direction  if  EOF  is  reached  first 


Entry  Points: 

None 

Elements: 

I/O 

Type 

Variable  Dimension 

Description 

I/O 

1 

IPR 

Profile  number 

1 

R 

Z 

Z value  of  last  point 

1 

R 

X 

X value  of  profile  to  store 

I 

1 

J 

Index  of  data  storage  location 

Comnxjn  Areas: 
Subroutines  Used: 
Link: 


/tlATA/,APlNT/,ATORE/,AlNlTS/,/CONST/,APDLPR/ 

disc,rwrite,pk 

STORE 


7.99  Subroutine  XMN  (X,YMX,YMN) 


Purpose: 
Entry  Points: 
E lements: 


Corrunon  Areas: 
Subroutines  Used: 
Link: 


To  solve  for  the  maximum  and  minimum  Y volue  given  an  X value 
None 

I/O  Type  Vorioble  Dimension  Description 

1 R X - X value  input 

0 R YMX  - Maximum  Y value 

0 R YMN  - Minimum  Y value 

Apdlpr/,Ajnits/ 

None 

PROFIL 


7.100  Subroutine  XY1NT  (X,YINTR,RMXY) 


Purpose: 


Computes  maximum  Y and  minimum  Y value  given  an  X . 
Also  computes  the  limits  in  X and  Y. 


Entry  Points: 
Elements: 


Common  Areas: 
Subroutines  Used: 
Link: 


None 


J/2 

tf 

Variable  Dimension 

Description 

1 

R 

X 

X value 

0 

R 

YINTR  2 

Y intersections  given  X value 

0 

R 

RMXY  4 

Maximum  and  minimum  X&Y  v 
for  particular  X value 

AH0T0/,/M1NMAX/,APDLPR/ 

MOVE 

CREATE 
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7.101  Subroutine  XYZPLH(XYZ,PLH) 


Purpose: 

To  convert  from  geocentric  coordinates  X,Y,Z  to  geographic 
coordinates  latitude,  longitude  (in  radians)  and  height  (in  meters) 

Entry  Points: 

None 

Elements: 

I/O  Type  Variable  Dimension  Description 

1 R XYZ  (3)  Geocentric  XYZ  in  meters 

0 R PLH  (3)  Geographic  latitude,  longitude 

ar>d  height 

Common  Areas: 

/GEODET/ 

Subroutines  Used: 

None 

Link: 

INPUT 

1 

Subroutine  YFILL  (Z,YI,l'»-X,l,IPR,X,YM,YlM) 

Purpose; 

To  fill  block  IF  maximum  Y value  on  tape  Is  encountered 
before  block  is  full 

Entry  Points: 

None 

Elements: 

I/O  Type  Variable  Dimension  Description 

1 R Z - Z value  of  last  recorded  point 

1 R Yl  - Y value  of  lost  recorded  point 

1 1 lYX  - Moximum  Y value  of  block 

III  - Location  of  last  stored  point 

In  common  array 

I/O  1 IPR  “ Profi  le  number 

1 R X - X value  of  profile 

I/O  R YM  - Minimum  Y value  to  store 

I/O  R YLM  - Maximum  Y value  to  store 

Common  Areas: 

/t)ATA/,APINT/,/STORE/,AJNITS/,/FPDLPR/,/CONST/ 

Subroutines  Used: 

PROPK,DISC,RWRITE 

Link: 

STORE 
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8. 


DESCRIPTION  OF  OOPS  INPUT  FILE 


The  OOPS  system  uses  three  terrain  data  sources  as  input  for 
the  generation  of  profiles  representing  the  terrain  of  the  rectified 
photograph.  The  contour  mode  uses  a MMS-32  data  tape  as  input.  The 
profile  mode  uses  a reformatted  AS-llB-1  profile  tape.  Input  for  the 
random  point  mode  is  a magnetic  tape  containing  point  ID  and  X,Y,Z 
ground  coordinates  of  discrete  points. 

8.1  Contour 

f1MS-32  data  tapes  contain  cultural,  hydrographic,  and  topo- 
graphic data.  These  different  data  types  are  identified  within  a header 
block.  The  x and  y coordinates  of  each  point  in  the  contour  data  tape 
are  contained  in  36  bit  words.  A complete  description  of  the  input 
information  for  the  contour  mode  can  be  found  in  the  Automatic  Carto- 
graphic System  Mod  III  by  Pennsylvania  Research  Associates,  Inc.  which 
was  developed  under  RADC  Contract  No.  RADC-TR-71-50. 

8. 2 Profile 

The  profile  mode  uses  as  input,  a data  tape  generated  on  the 
PDP  11/45.  This  tape  is  the  result  of  processing  an  AS-llB-1  data  tape 
through  a program  designed  to  take  out  any  irregularities  in  the  data. 
The  magnetic  tape  is  blocked  in  256,  16  bit  words.  Generation  of  this 
tape  is  performed  by  running  a program  called  GRIDIN.  This  program 
was  written  by  RADC  personnel. 
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The  magnetic  tape  generated  from  the  GRIDIN  program  is 
blocked  in  256  integer  words  or  128  real  words.  All  x,y  and  z 
coordinates  are  floating  point  numbers.  The  format  of  each  block 
is  as  follows:  real  word  1 and  2 is  blank,  words  3-44  contain  x 

values,  words  45-86  contain  y values  and  words  87-128  are  the  z 
values  associated  with  the  x and  y coordinates.  The  x coordinates 
are  incremented  starting  at  the  minimum  x value  in  1 millimeter 
increments.  The  y value  is  incremented  every  250  microns.  All 
coordinates  are  given  in  millimeters. 

8. 3 Random  Points 

Input  to  the  random  point  mode  consists  of  a magnetic  tape 
containing  the  Identification  number  and  ground  coordinates  of  dis- 
crete points.  This  data  is  written  with  an  unformatted  write  con- 
sisting of  one  integer  word  (ID)  and  three  real  words  (X,Y,Z).  Each 
point  is  written  in  this  manner.  The  file  is  ended  when  a point  is 
input  with  an  identification  number  of  9999. 
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9. 


DESCRIPTION  OF  OOPS  OUTPUT  FILE 


The  generation  of  a rectified  orthophotograph  on  the  offline 
orthophoto  printer  is  accomplished  by  using  two  input  files.  The  first 
of  these  files  is  the  decimal  shutdown  tape.  This  paper  tape  contains 
the  information  necessary  for  orientation  of  the  photograph  on  the  OOP. 
The  second  file  is  a magnetic  tape  containing  the  data  representing  the 
terrain  of  the  rectified  orthophotograph.  Both  tapes  duplicate  the 
format  of  data  output  from  the  AS-llB-1. 

9. 1 Decimal  Shutdown  Tape 

The  decimal  shutdown  tape  contains  all  the  constants  necssary 
for  the  orientation  of  the  photograph  on  the  OOP  along  with  an  alpha- 
numeric identifier  that  precedes  each  data  constant.  Thus  shutdown 
tape,  in  the  form  of  a paper  tape,  is  punched  in  RCA  501  code.  All 
information  punched  on  this  paper  tape  is  also  tabulated  on  hard  copy. 

Information  pertaining  to  orientation  is  in  an  AS-11  coordin- 
ate frame.  This  coordinate  frame  uses  the  projected  principle  point  as 
the  origin  and  the  +Y  photo  axis  as  the  +Y  model  axis  (see  Section  4.1). 
The  data  constants  contained  on  the  shutdown  tape  include  values  for 
both  fore  and  aft  photographs.  This  allows  the  user  to  mount  either 
photo  for  rectification.  Data  constants  include  the  offsets  in  x and 
y,  focal  length,  radius  of  the  earth,  point  of  tangency  of  the  photo, 
and  constants  for  atmopsheric  refraction.  Nadir  point  values,  air  base 


components  and  linear  rates,  and  orientation  and  orientation  rates  are 
also  included.  In  addition  photograph  type  (pan  or  frame)  is  available 
on  this  shutdown  tape. 

9.2  OOPS  Profile  Tape 

The  OOPS  profile  tape  contains  the  derived  terrain  information. 
This  tape  is  formatted  in  18-bit  integer  words  of  "X,Y,Z"  coordinates 
and  an  associated  identification  label.  Between  any  two  variables, 
there  is  an  18-bit  word  with  a value  of  zero.  These  18-bit  integer 
words  are  derived  from  one  POP  real  word.  The  conversion  from  16  to 
18-bit  words  is  accomplished  by  computing  the  binary  representation  of 
the  real  POP  word.  This  binary  representation  is  then  output  to  a 
seven  track  magnetic  tape  to  be  input  on  the  BX-272. 

The  output  tape  is  compatible  with  the  input  requirements  of 
the  off-line  orthophotoprinter.  Following  is  a diagram  of  the  output 
tape  format: 


3S" 

3its 

W 

18  bits |18  bits 

18  bits 

18  bits 

18  bits 

18  bits 

BX-272  wds 

Integer 

Integer 

Interger 

Integer 

Word  Type 

Word  Content 

Zero  1 ID 

atm 

Zero 

Y 

Zero 

Z 
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The  tape  is  blocked  in  physical  records  of  256  BX-272  18-bit 
words.  Thus,  32  data  points  are  contained  in  each  physical  record.  The 
initial  profile  is  increasing  in  the  "Y"  coordinate  while  the  second 
profile  decreases  in  "Y".  The  data  is  compacted  between  profiles  such 
that  data  from  two  profiles  can  be  contained  in  one  physical  record. 
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10. 


RESULTS 


Considerable  testing  using  both  real  and  simulated  data  was 
accomplished  during  the  OOPS  contract.  However,  time  did  not  permit 
extensive  compilation  of  statistics.  Both  real  and  simulated  tests 
were  performed  for  the  three  program  modes  (contour,  profile  and  ran- 
dom points)  and  for  the  coordinate  transformations  needed  to  generate 
the  shut-down  tape.  The  basic  approach  was  to  compare  plots  of  the 
real  data  inputs  with  plots  of  the  generated  OOPS  output. 

The  simulated  data  sets  were  all  easily  recognizable  geometric 
figures  such  that  any  deviation  from  that  geometric  shape  was  readily 
discernible.  For  the  contour  and  profile  modes,  two  hemispheres  placed 
side  by  side  were  used  as  fictitious  data.  These  hemispheres  were  of 
different  heights,  with  one  hemisphere  ranging  from  a height  of  1,000 
feet  to  2,500  feet  and  the  second  from  1,000  feet  to  1,500  feet.  For 
the  random  points  mode,  a surface  was  created  from  a sine  wave  function. 
Each  profile  of  this  surface  was  an  elongated  sine  curve  with  each 
successive  profile  being  translated  parallel  to  the  previous  one. 

The  real  data  tests  were  all  developed  over  the  South  Mountain 
area  near  Phoenix,  Arizona.  This  well  known  test  area  is  noted  for  its 
sharp  contrast  in  terrain  relief.  Slope  in  this  area  varies  from  ap-  . 
proximately  0°  to  in  excess  of  60®,  which  presents  a difficult  problem 
of  terrain  modeling.  The  following  paragraphs  describe  the  results 

obtained  from  each  of  the  three  modes  and  the  coordinate  transformations. 
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ContouA  Modz.  Figures  10.1-10.5  illustrate  the  simulated 


data  inputs,  results  from  simulated  data,  actual  data  inputs,  and  re- 
sults from  real  data  inputs,  respectively  for  the  contour  mode.  The 
simulated  hemispheres,  as  shown  in  Figure  10.1,  were  digitized  in  MMS- 
32  format  by  RADC.  This  digitized  data  was  converted  to  profile  infor- 
mation and  the  results,  as  output  from  the  OOPS  software,  are  illustra- 
ted in  Figure  10.2.  Excellent  results  were  obtained  from  the  contour 
mode  over  areas  of  moderately  sloping  relief  using  a polynomial  func- 
tion, however,  an  interpolation  routine  provided  more  successful  results 
over  extremely  rugged  terrain.  Figure  10.4  illustrates  results  obtain- 
ed from  the  real  data  input  set  in  Figure  10.3  using  a polynomial  func- 
tion to  derive  the  profiles.  These  results  are  less  accurate  due  to 
the  sparsely  digitized  contours  and  extreme  magnitude  of  the  terrain 
slope  changes.  However,  when  compared  with  the  actual  AS-llB-1  profiles , 
shown  in  Figure  10.7,  these  profiles  do  conform  to  the  actual  terrain 
in  most  areas.  The  same  input  data  was  used  to  generate  profiles  using 
an  interpolation  routine.  As  illustrated  in  Figure  10.5,  this  method 
of  generating  profiles  from  rugged  areas  more  closely  resemble  the  actual 
terrain.  Figures  10.4  and  10.5  were  plotted  at  a different  scale  and 
orientation  than  the  actual  AS-llB-1  profiles  shown  in  Figure  10.7.  Gen- 
eral terrain  characteristics  can  be  compared,  however,  a direct  compari- 
son between  the  three  figures  cannot  be  made. 
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I 

Mode.  The  profile  mode  was  tested  in  a very  similar 
manner.  The  profiles,  as  derived  from  the  hemispheres  in  the  contour 
mode,  were  used  as  AS-llB-1  input  data  to  the  profile  mode.  This  simu- 
lated data  was  processed  through  the  profile  mode  and  the  results  are 
illustrated  in  Figure  10.6.  This  plot  of  a small  portion  of  the  sphere 
illustrates  an  inherent  problem  in  modeling  smooth  surfaces  with  the 
multiquadric  function.  Each  of  the  small  bumps  on  the  spherical  sur- 
face are  nodal  points  about  which  the  multiquadric  function  was  deri- 
ved. Recall,  Figure  10.7  is  a plot  of  the  original  profile  data  as 
derived  manually  from  the  AS-llB-1.  The  profiles  derived  by  the  OOPS 
software  were  oriented  to  correspond  to  that  of  the  input,  so  that  a 
comparison  could  be  made,  however,  these  output  profiles  can  be  oriented 
to  any  desired  azimuth  angle.  The  output  from  the  profile  mode  is  illus- 
trated in  Figure  10.8.  These  results  are  felt  to  be  adequate  for  ortho- 
photo production.  Some  features  are  generalized  due  to  the  multiquadric 
functional  representation.  This,  however,  is  quite  desirable  for  ortho- 
photo compilation. 

Random  Po^nti.  Two  data  sets  were  processed  using  the  random 
points  mode.  The  first  was  a simulated  data  case  where  scattered  points 
were  generated  from  a three-dimensional  sine  wave  function.  The  second 
data  set  was  real  data  which  consisted  of  discrete  points  measured  from 
a map  sheet  in  a UTM  coordinate  frame.  The  real  data  set  again  covered 
the  South  Mountain  area.  Plots  are  not  available  of  the  generated 
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profiles  for  either  the  simulated  or  real  data  cases.  However,  these 
profiles  were  generated  from  the  same  mathematical  expression  as  that 
used  in  the  profile  mode.  The  results  from  the  sine  wave  function  were 
consistent  and  seemed  to  approximate  the  original  sine  function  quite 
well.  The  resultant  RMS  value  for  all  points  in  the  simulated  area  was 
18  meters.  Over  the  real  data  area  262  data  points  were  available, 
with  an  RMS  value  for  these  points  of  35  meters. 

CooA.cU.nate.  T^n&^omcvtioni.  Coordinate  transformations  de- 
veloped during  this  contract  for  the  generation  of  the  footprint  boun- 
daries and  shut-down  tape  were  validated  using  real  data  from  previous 
DBA  frame  and  panoramic  reductions.  Transformations  were  developed  re- 
lating USR,  LSR,  LV,  Geographies  and  UTM  coordinates.  Photograph  ori- 
entations which  were  converted  from  one  coordinate  system  to  another 
were  verified  by  projecting  known  image  points,  using  the  transformed 
orientation  matrix,  and  comparing  the  resultant  object  space  coordinates 
to  their  known  values. 

Routines  were  developed  for  computing  footprint  boundaries  or 
ground  coordinates,  by  projecting  image  coordinates  to  the  earth.  These 
routines  were  tested  using  both  frame  and  panoramic  image  coordinates 
with  known  ground  coordinates. 
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FIGURE  10.3.  Contour  Input,  South  Mountain  Area 


FIGURE  10.7.  AS-llB-1  Original  Profile  Data,  South  Mountain  Area 


11. 


CONCLUSIONS  AND  RECOMMENDATIONS 


Several  methods  of  modeling  terrain  relief  were  investigated 
and  implemented  during  the  OOPS  contract.  These  methods  included; 
multi  quadric,  polynomial,  interpolation  and  a combination  of  inter- 
polation and  multiquadric.  For  patterned  input  data,  such  as  the  pro- 
file data,  the  multi  quadric  function  proved  to  be  a satisfactory  method 
of  representing  the  terrain.  For  unpatterned  data,  three  of  the  above 
mathematical  approaches  were  developed.  These  were:  polynomial,  inter- 

polation and  a combination  of  interpolation  and  multiquadric.  Random 
points  are  processed  using  a multiquadric  function,  which  provided  to 
be  quite  adequate.  Contour  data  is  processed  using  both  a polynomial 
function  and  an  interpolation  routine.  The  results  achieved  using  the 
polynomial  function  are  excellent  for  areas  of  moderately  sloping  ter- 
rain, however,  in  areas  of  steep  terrain,  an  interpolation  routine 
provided  the  best  results.  These  mathematical  approaches  are  applied 
to  data  of  various  types  to  provide  data  in  a format  suitable  for  the 
off-line  orthophoto  printer.  Orientation  of  the  final  model  can  be 
derived  from  nearly  any  input  coordinate  system,  including  USR,  LSR, 

LV,  Geographies  and  AS-11  absolute  model  parameters.  Data  used  for 
the  generation  of  output  profiles  can  be  computed  from  AS-llB-1  pro- 
files, MMS-32  contour  data  and  random  points.  Combinations  of  the 
above  data  types  are  also  acceptable.  This  comprehensive  software 
system  provides  a major  contribution  for  expanding  the  capability  of 
the  off-line  orthophoto  system. 
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The  OOPS  system  was  designed  and  implemented  primarily  as  a 
research  and  development  program.  Several  modifications  to  this  system 
could  be  made  to  improve  flexibility  and  computation  efficiency.  Em- 
ploying a combination  of  array  algebra  and  the  multiquadric  function 
in  the  profile  and  random  point  modes  would  improve  software  efficiency. 
Considering  the  unpatterned  data  distribution  in  the  random  point  mode, 
an  interpolation  routine  would  be  required  prior  to  beginning  the  array 
algebra  solution.  This  interpolation  routine  would  establish  grid 
points  having  x and  y coordinates  and  an  interpolated  height.  The  in- 
terpolated height  of  any  grid  point  would  be  computed  from  a function 
which  utilizes  distances  and  heights  of  surrounding  data  points.  Be- 
cause the  input  profile  data  is  already  in  a gridded  format,  no  data 
reformatting  is  required  for  that  mode.  Array  algebra  could  greatly 
reduce  both  the  time  and  computer  core  required  to  compute  the  coef- 
ficients representing  terrain  relief.  At  the  present  time,  the  solu- 
tion of  an  area  70mm  x 110mm  with  a 6x6  nodal  pattern  requires  approxi- 
mately 3.5  hours.  Modification  of  the  software  to  model  the  multi - 
quadric  function  using  array  algebra  could  easily  reduce  this  computa- 
tional time  by  a factor  of  ten.  It  is  strongly  recommended  that  the 
OOPS  system  be  modified  so  that  the  solution  for  the  multiquadric 
function  be  performed  using  array  algebra. 

A second  improvement  to  OOPS  would  be  to  determine  and  use 
the  "best"  nodal  pattern  and  constant  value  for  each  different  type  of 

terrain,  whenever  the  multiquadric  math  representation  is  employed. 
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The  determination  of  these  two  items  could  be  performed  either  manually 
or  automatically.  This  improvement  would  result  in  time  savings  as 
well  as  increased  accuracy.  Areas  with  fairly  level  terrain  require  a 
sparce  nodal  pattern;  the  number  of  computations  become  less  and  the 
size  of  the  matrices  in  core  are  smaller.  When  comparing  a nodal  point 
pattern  of  3x3  with  one  of  7x7,  a significant  amount  of  computation 
time  is  saved.  There  are  approximately  one-fifth  as  many  computations 
in  a 3x3  nodal  pattern  as  compared  to  a 7x7  nodal  pattern.  Thus,  one 
could  select  the  optimum  nodal  pattern  to  achieve  maximum  efficiency 
while  maintaining  adequate  accuracy. 

The  modular  design  of  the  OOPS  system  could  be  expanded  to 
accept  new  data  types.  One  data  type  which  presently  exists  and  was 
not  addressed  during  this  contract  is  LIS.  This  data  consists  of  vec- 
tor rays  between  data  points.  Other  data  types  which  could  be  in- 
corporated into  the  OOPS  system  are  TERCOM  and  DRLMS.  Such  data  is 
given  in  an  evenly  gridded  format  and  therefore  would  be  ideally  suit- 
ed to  the  multiquadric  function  solution  using  array  algebra.  The 
ability  to  create  profiles  using  any  of  the  above  data  types  would 
greatly  increase  the  capability  of  the  OOP.  Another  improvement 
which  could  be  incorporated  in  the  OOPS  is  the  ability  to  read  an 
AS-llB-1  data  tape  created  directly  from  the  stereo  plotter.  Present- 
ly this  raw  data  is  input  to  a separate  software  system  for  editing 
prior  to  entry  into  OOPS.  This  editing  function  could  be  accomplished 

directly  in  OOPS  without  an  intermediate  software  step. 
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The  major  obstacle  encountered  during  the  OOPS  development 
was  the  extremely  large  computer  storage  requirement.  The  core  re- 
quirement resulted  from  both  the  large  amount  of  input  data  to  be 
processed  and  the  software  code  itself.  The  available  mass  storage 
was  also  very  limited  since  only  one  RK05  disc  is  available  on  the 
RADC  POP  11/45.  This  created  several  I/O  problems  and  considerable 
time  was  lost  in  cleaning  and  compacting  disc  files.  It  is  recommen- 
ded that  a minimum  of  one  more  RK05  disc  be  purchased.  During  this 
contract,  files  had  to  be  limited  and  therefore  data  sizes  were  also 
limited.  Much  of  the  I/O  currently  done  on  magnetic  tape  could  be 
accomplished  more  efficiently  on  disc.  Further  recommendations  as  to 
hardware  would  include  the  purchase  of  the  RSX-llM  operating  system 
with  Fortran  IV  Plus.  The  purchase  of  RSX-llM  would  allow  the  user 
to  access  the  full  64K  of  core.  Fortran  IV  plus  would  enable  the 
programmer  to  use  such  features  as  subroutine  entry  points  and  double 
word  integers. 

During  any  potential  upgrade  of  the  OOPS,  it  is  recommended 

that  the  individual  processes  be  broken  up  into  subtasks  under  control 

of  a RSX  monitor.  Between  these  subtasks,  checkpointing  or  operator 

interaction  is  recommended.  Here  the  operator  can  analyze  the  output 

at  such  points  and  either  rerun  that  part  of  the  program  or  continue. 

This  modification  would  result  in  a more  efficient  software  system  in 

terms  of  utilization  of  computer  operation  time.  After  debugging  a 

complete  mode,  such  as  contour,  profile  or  random  points  reduction, 
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the  programs  should  be  tested  using  several  data  cases  to  exercise 
each  option.  After  all  modes  have  been  debugged  in  this  manner,  the 
complete  system  could  be  incorporated.  Each  system  link  of  the  OOPS 
software  requires  approximately  35  minutes  of  run  time.  Therefore, 
much  time  could  be  lost  by  not  debugging  programs  as  modules.  How- 
ever, although  development  and  "debug"  is  easier  when  the  software 
is  in  separate  modules,  eventual  utilization  is  tedious  unless  the 
individual  modules  are  combined  under  a single  controller. 

In  summary,  the  development  of  the  OOPS  has  enhanced  the  cap- 
ability of  the  off-line  orthophoto  printer  in  that  various  types  of 
data  can  now  be  used  for  the  computation  of  profile  tapes  for  the  Off- 
line Orthophoto  Printer.  The  current  software  system  demonstrates  the 
feasibility  of  the  multiquadric  and  polynomial  functions  to  represent 
terrain  surfaces  with  acceptable  accuracies.  Incorporation  of  the 
above  recommendations  into  the  OOPS  would  result  in  a substantial 
improvement  in  run  time,  as  well  as  an  increased  product  accuracy. 
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